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REVISION HISTORY 


This manual is the APMATH64 Manual, Volume 3, 869-7482-991. The letter 
shown under the revision number column indicates the portion of the 
part number that changes for each revision. The last entry is the 
latest revision to this manual. 


The revision history begins with this manual. 


Deleted Utilities Library, deleted the 
LPSPFI subroutine, added internal subroutine 


information, and added 16 new routines. 1/87 

Added routines to Basic Math Library 

Double Precison Library, and Matrix 

Algebra Accelerated Math Library. 12/87 
NOTE: For revised manuals, a vertical line "|" outside the left 


Margin of the text signifies where changes have been made. 


NOTE TO READER 


This is the third volume of the APMATH64 Manual. It 
is comprised of part 3 of Appendix A and Appendix B 
through Appendix J. Note that Appendix A. continues 
through Volumes 1, 2, and 3. The page numbers are 
listed consecutively through the volumes. 


The APMATH64 Manual has three indices located at the 
end of Volume 3 and two at the end of Volume 4. The 
first index (Appendix I) is a list of the APMATH64 
routines in page order by type. The second index 
(Appendix J) is an alphabetical list of all the 
APMATH64 routines. The third index is a key word 
index of the APMATH64 routines. The fourth index 
(Appendix L) is an alphabetical list of the 
APMATH64/MAX routines. The fifth index is a key word 
index of the APMATH64/MAX routines. 
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EREEREEKER RAEREEREERE 
= * & & 
* MMTMUL * ——— VECTOR MULTIPLY (MD*MD TO TM) —-- * MMTMUL * 
& * x & 
KKEKEKEERE PEREEREAELE 
PURPOSE: To multiply the elements of two vectors in Main 

Memory and store the resultant vector in Table 

Memory. 


CALL FORMAT: CALL 


PARAMETERS : 


5 


i 
" 
‘ 
¢ 


SAY owoH PY 


MMTMUL(A,I,B,J,ITMC,K,N) 


= Floating-point Main Memory input vector 
= Integer element step for A 
= Floating-point Main Memory input vector 


= Integer element step for B 
= Integer base address of TM output vector C 


ess eeyé Seer ae SS 


= Integer element step for C 
= Integer element count 


DESCRIPTION: MMTMUL multiplies N elements of the vector A with N 
elements of the vector B, where A and B are in Main 
Memory, and stores the results in a vector with base 
address ITMC and increment K in Table Memory. 


NOTE: Writable Table Memory begins at address 8192. 


EXAMPLE: 
N=3 
I=J=K=1 
ITMC = 8192 
A : 1.9 2.0 3.9 
3 3.8 4.2 5. 


TMLOC: 8192 8193 8194 
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: 3.0 8.8 15.9 
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REREKREERE RREREEEEEE 
x e 3 & & 
* MTIMOV * ——— VECTOR MOVE WITH INCREMENT (MD TO TM) ——- * MTIMOV * 
* = * & 
RREKKREKEER RPERERRARERE 
PURPOSE: To move elements of a vector from Main Memory 


to Table Memory, where the increments between the 
elements are specified. 


CALL FORMAT: CALL MTIMOV(A,I,ITMC,K,N) 
PARAMETERS : = Floating-point Main Memory input vector 

= Integer element step for A 
TMC = Integer base address of TM output vector C 


= Integer element step for C 
= Intecer element count 


aoe eS —— > = See ee eS fe 


AA AH Pp 


DESCRIPTION: MTIMOV moves the elements of an input vector A with 
increment I in Main Memory to an output vector with 
base address ITMC and increment K in Table Memory. 


NOTE: Writable Table Memory begins at address 8192. 


EXAMPLE: 
N = 3 
K = 2 
ITMC = 8192 


A >: 1.8 2.9 3.9 


TMLOC: 8192 8193 8194 8195 8196 8197 
Cc >: 1.8 X 2.8 X 3.9 X 


X represents unchanged values. 
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2 & & & 
* MIMMUL * ——-— VECTOR MULTIPLY (MD*TM TO MD) -—- * MIMMUL * 
; J 7 2 * 
RENREKAEARAEEE RERAEREAETTE 
PURPOSE: To multiply elements of a vector in Main Memory 


by elements of a vector in Table Memory and 
store the products in Main Memory. 


CALL FORMAT: CALL MTMMUL(A,1I,ITMB,J,C,K,N) 


PARAMETERS : 


TMB = 


SAO UHH YP 
\ 


DESCRIPTION: MTMMUL 
Memory 


Floating-point Main Memory input vector 
Integer element step for A 
Integer base address of TM input vector B 


= Integer element step for B 


Floating-point Main Memory output vector 


Integer element step for C 
Integer element count 


multiplies N elements of the vector A in Main 
by N elements of the vector with base address 


ITMB in Table Memory, and stores the products in N 
@élements of the vector C in Main Memory. 


EXAMPLE: 
N=3 
I=J=K=1 
ITMB=8192 
A : 1.9 2.0 3.8 
TMLOC: 8192 8193 8194 
B : 2.0 3.2 4.9 
Cc H 2.9 6.8 12.9 
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RREREUEAREE RREAEEEREEK 
& = & * 
* MTMSUB * -——- VECTOR SUBTRACT (MD-TM TO MD) —~ * MTMSUB * 
& . 4 . x 
RRKEKREEEE RREREEKEERE 
PURPOSE: - To subtract the elements of a vector in Table 


Memory from the elements of a vector in Main 
Memory and store the results in a vector in 
Main Memory. 


CALL FORMAT: CALL MTMSUB(A,I,ITMB,J,C,K,N) 


PARAMETERS : A = Floating-point Main Memory input vector 
L = Integer element step for A 
ITMB = Integer base address of TM input vector B 


J = Integer element step for B 
¢ = Floating-point Main Memory output vector 
K = Integer element step for C 
N 


= Integer element count 
DESCRIPTION: MTMSUB subtracts N elements of a vector with base 
address ITMB in Table Memory from N elements of the 


vector A in Main Memory, and stores the results in WN 
elements of the vector C in Main Memory. 


EXAMPLE: 


ITMB = 8192 
A : 3.9 4.9 5.8 


TMLOC: 8192 8193 8194 
B : 2.0 1.2 1.9 


C : 1.0 3.9 4.9 
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PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 
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RERKEEKEEE 


2 * 


* MPTMUL * 
& & 


-—— VECTOR MULTIPLY (MD*TM TO TM) —— 


REREEEEKAEK 


To multiply the elements of a vector in Main 
Memory by the elements of a vector in Table 
Memory and store the products in a vector in 
Table Memory. 


CALL MTTMUL(A,I,ITMB,J,ITMC,K,N) 


A = Floating-point Main Memory input vector 

a = Integer element step for A 

ITMB = Integer base address of TM input vector B 
J = Integer element step for B 

ITMC = Integer base address of TM output vector C 
K = Integer element step for C 

N = Integer element count 


MTTMUL multiplies N elements of the vector A in Main 
Memory by N elements of the vector with base address 
ITMB in Table Memory, and stores the products in N 
elements of a vector with base address ITMC in Table 
Memory. 


NOTE: Writable Table Memory begins at address 8192. 


N=3 

I=J=K=1 

ITMB = 8192 

ITMC = 8292 

A : 3.8 4.9 5.2 
TMLOC: 8191 8193 8194 
B 2.9 1.9 3.8 
TMLOC: 8292 8293 8294 
C : 6.9 4.9 15.9 
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* TMDOT * 
* *® 
RRRAEEREKE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 
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——- REAL DOT-PRODUCT (TM AND MD) -——— 


Computes the real dot-product of two vectors 
where one vector is stored in Main Memory and 
the other vector is stored in Table Memory. 


APPENDIX A 


RURAEKAKRERE 
= & 
* T™DOT * 
= . 
REREEEEREE 


Both vectors are assumed to be stored compactly. 


CALL TMDOT (ITMA,B,C,N) 


ITMA = Integer base address of TM input vector A 


= Floating-point Main Memory input vector 
Floating-point Main Memory output scalar 


20 W 
i) 


TMDOT computes the real dot-product of N elements of 
the vector with base address ITMA in Table Memory with 
N elements of the vector B in Main Memory, and stores 


the resultant scalar in Main Memory. 


Formula: 
C = A(1)*B(1) + A(2)*B(2) + ... + A(N)*B(N) 
C= 8.8, if N<l 


N = 3 
ITMA = 8192 
TMLOC 8192 8193 8194 
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* MMM * 
, x 
REEREKEEETE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 


APPENDIX A 


RekKRERRES 
* * 
-——~- MATRIX MULTIPLY (TM WORKSPACE) --—- * TMMM * 
* * 
REKRARAEEE 


Multiplies two matrices A and B in Main Memory 
to form a matrix C in Main Memory. This routine 
uses a workspace in Table Memory to achieve high 
speed. 


CALL TMMM (A,B,C,MC,NC,NA, ITMW) 


A = Floating-point Main Memory input matrix 
B = Floating-point Main Memory input matrix 
Cc = Floating-point Main Memory output matrix 
MC = Integer number of rows in output matrix C 


(and input matrix A) 


NC = Integer number of columns in output matrix C 
(and input matrix B) 
NA = Integer number of columns in input matrix A 


(and number of rows of input matrix B) 
ITMW = Integer base address of TM work area of length NA 


TMMM computes the product of the MC-row by NA-column 
matrix A and the NA-row by NC-column matrix B (both in 
Main Memory) and stores the result in the MC-row by 
NC-column matrix B in Main Memory. This routine 

uses a workspace of length NA in Table Memory to achieve 
high speed. All matrices are assumed to be stored in 
column order. 


NOTE: Writable Table Memory begins at location 8192. 
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RERRERREEE REKKEEKEERE 
& . & = 
* TMMSUB * -——— VECTOR SUBTRACT (TM-MD TO MD) ——— * TMMSUB * 
x & z * 
REKREEREREE REKKEKRERE 
PURPOSE: To subtract the elements of a vector in Main 


Memory from the elements of a vector in Table 
Memory and store the results in a vector in 
Main Memory. 


CALL FORMAT: CALL TMMSUB(ITMA,1I,B,J,C,K,N) 


PARAMETERS : ITMA = Integer base address of TM input vector A 
= Integer element step for A 

= Floating-point Main Memory input vector 

= Integer element step for B 

Floating-point Main Memory output vector 
= Integer element step for C 


= Integer element count 


ZRaaoaaz 
rT) 


DESCRIPTION: TMMSUB subtracts N elements of the vector B in Main 
Memory from N elements of the vector with base address 
ITMA in Table Memory, and stores the differences in the 
vector C in Main Memory. 


EXAMPLE: 
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® & & & 
* TMVLC2 * ~—— VECTOR LINEAR COMBINATION --—— * TMVLC2 * 
* e 2 5 & 
RERXEREREER REKEREKATRE 
PURPOSE: To compute the linear combination of two vectors, 


one in Table Memory and the other in main memory, 
and store the resultant vector in main memory. 


CALL FORMAT: CALL TMVLC2 (Sl, ITMA, S2, B, J, C, K, N) 


PARAMETERS : Sl = Floating-point scalar coefficient 
for the TM input vector A 
ITMA = Integer base address of the TM input 
vector A 


e9 = FI es - > Teh 
S2 = Floating-point scalar coefficient 


for the MD input vector B 
B = Floating-point MD input vector 
J = Integer element step for B 
Cc = Floating-point MD output vector 
K = Integer element step for C 
N = Integer element count 


DESCRIPTION: C(m) = Sl * A(m) + S2 * B(m); for m= 1toN 


Where A is in Table Memory, and B, Sl, S2, and C 
are in main memory. 


EXAMPLE: 

N = 3 

Sl = -1.9 

S2 2.8 

J = 1 

K = 1 

ITMA = 8192 

TMLOC: 8192 8193 8194 
A : 1.8 2.8 3.8 
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REREEEEEER 3 REREEREKEE 
. J x & * 
* TTMADD * —— VECTOR ADD (TM+TM TO MD) -——— * TTMADD * 
& . 7 * 
RRERETZTEKRR REeKERXERETE 
PURPOSE: To add the elements of two vectors in Table 


Memory and store the sums in Main Memory. 
CALL FORMAT: CALL TTMADD(ITMA,1I,ITMB,J,C,K,N) 


PARAMETERS : ITMA = Integer base address of TM input vector A 
r = Integer element step for A 
ITMB = Integer base address of TM input vector B 
J = Integer element step for B 
Cc = Floating-point Main Memory output vector 
K 
N 


to 410 YCL oo we oth eS i BMS dl LVL ed 


= Integer element count 


DESCRIPTION: TTMADD adds N elements of the vector with base address 
ITMA in Table Memory to N elements of the vector with 
base address ITMB in Table Memory, and stores the sums 
in N elements of the vector C in Main Memory. 


EXAMPLE: 
N=3 
T=J=K=1 
ITMA = 8192 
ITMB = 8292 


‘) 


5.2 7.6 9.6 
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RREKEEREERE RERETKEEEEE 
x & x x 
* TTMSUB * ——— VECTOR SUBTRACT (TM—-TM TO MD) —— * TTMSUB * 
; 2 > : > & 
RRRRRERREE ERREKRARERE 
PURPOSE: To subtract the elements of two vectors in Table 


Memory and store the differences in a vector in 
Main Memory. 


CALL FORMAT: CALL TITMSUB(ITMA,I,ITMB,J,C,K,N) 


PARAMETERS : ITMA = Integer base address of TM input vector A 
I = Integer element step for A 
ITMB = Integer base address of TM input vector B 
= Integer element step for B 
= Ploatina-coint Main Memory output vector 
= Integer element step for C 
= Integer element count 


DESCRIPTION: TTMSUB subtracts N elements of the vector with base 
address ITMB in Table Memory from N elements of the 
vector with base address ITMA in Table Memory, and 
stores the resulting differences in a vector C in Main 
Memory. 


EXAMPLE: 
N=3 
I=J=K=1 
ITMA = 8192 
ITMB = 8292 


TMLOC: 8192 8193 8194 


a 
re 
as] 
Ww 
M 
> 
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* TTTMUL * 
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REEXKEEREEE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 


—-— VECTOR MULTIPLY (TM*TM TO TM) —— 


APPENDIX A 


RREEEAERELX 
7 x 
* TTTMUL * 
* x 
REERVTAEEREE 


To multiply the elements of two vectors in Table 


vector in Table Memory. 


resulting products in a 


CALL TTTMUL(ITMA,I,ITMB,J,ITMC,K,N) 


TTTMUL 


NOTE: 


N=3 


I=J=K=1 


ITMA = 
ITMB = 
ITMC = 


TMLOC : 
A : 
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= Integer 
= Integer 


Integer 
Integer 


me am 


= Tak 
= integer 


Integer 
Integer 


base address 
element step 
base address 
element step 
| _ —_ mm 


| 
Mast auuLlLess 


element step 


element count 


of TM 
for A 
of TM 
for B 


me rine 
Wh Li 


for C 


input vector A 


input vector B 


wan 


multiplies N elements of the vector with base 
address ITMA in Table Memory by N elements of the 
vector with base address ITMB in Table Memory, and 
stores the resultant products in the vector with base 
address ITMC in Table Memory. 


8192 
8292 
8392 


8192 
1.8 


8193 8194 
2.0 3.8 
8293 9294 
4.0 5.9 
8393 3394 
8.0 15.9 


Writable Table Memory begins at address 8192. 
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x ed & * 
* TTVLC2 * -—-—— VECTOR LINEAR COMBINATION ——- * TTVLC2 * 
x & & x 
REKEKERKERE KREEEERERELR 
PURPOSE: To compute the linear combination of two vectors, 


one in Table Memory and the other in main memory, 


Wate ane 


and store the resultant vector in Table Memory. 
CALL FORMAT: CALL TTVLC2 (Sl, ITMA, S2, B, J, ITMC, N) 
PARAMETERS : Sl = Floating-point scalar coefficient for the 


TM input vector A 
ITMA = Integer base address of the TM input 


vector A 

SZ = Floating-point scalar coefficient for the 
MD input vector B 

B = Floating-point MD input vector 

J = Integer element step for B 

ITMC = Integer base address of the TM output 
vector C 

N = Integer element count 

DESCRIPTION: C(m) = Sl * A(m) + S2 * B(m); for m= 1toN 


Where A and C are in Table Memory, and B, Sl, 
and S2 are in main memory. 


Note: Writable Table Memory begins at address 8192. 


EXAMPLE: 
N 3 
Sl = -1.9 
S2 = 2.9 
J = 1 


TMLOC: 8195 8196 8197 
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RERERERERE 
ci = 


* PEEK * 
= 2 
RAREAREEERE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 


APPENDIX A 


RREKRAEKEKE 
PJ e 
——— MEMORY FETCH —— * PEEK * 
2 = 
ReEKKEEEKESE 


To fetch the contents of a specified memory word. 
Function Value = PEEK(Addr) 


Function Value = The unformatted contents of the 
specified memory location 

Addr = An integer specifying the address 
to be accessed 


The specified memory location is accessed and its 
contents returned as the function-valiue output. The 
output is the unformatted word. That is, no format 
conversion is performed by the function. 


(Assuming location 18@@ contains 
91 23 34 56 78 9A BC DE (hex) ) 


Addr 
Function Value 


1929 
G1 23 34 56 78 9A BC DE 
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x * & *& 

* VIFIX * ——~- VECTOR INTEGER FIX ——- * VIFIX * 

* * ®& & 

RERERRIEKE SREAKKERERE 

PURPOSE: To fix to 53-bit integers the elements of a 
floating-point vector. 


CALL FORMAT: CALL VIFIX(A,I,C,K,N,F) 
PARAMETERS: = Floating-point input vector 

= Integer element step for A 

= Long-integer output vector 

Integer element step for C 

= Integer element count 

= Integer flag (@ to round, i to truncate) 


myaROH YP 
iT) 


DESCRIPTION: C(m)=FIX(A(m)); for m=1 to N 


EXAMPLE: 
N=4 
F= g 
AM: 227 =1.5° =332. 3:5 
Coe. 22 —2 -3 4.9 
N= 4 
Fe=l 
A: 1.7 -1.5 -3.2 3.5 
C. # -1 -3 3.2 
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2 . x x 
* YPK16 * -——— VECTOR 16-BIT BYTE PACK -—— * YVPK16 * 
x . 3 & & 
RRRRREEEEE RRREEKKERE 
PURPOSE: To pack each four 64-bit floating-point numbers 

into one destination word as 16-bit quarter words 


CALL FORMAT: CALL VPK16(A,I,C,K,N,F) 


PARAMETERS : Floating-point input vector 

= Integer element step for A 

= Signed-quarterword~integer output vector 

= Integer element step for C 

= Integer element count (destination words) 
) 


Integer flag (2 to round, 1 to truncate 


ym ROH YP 
Q 


DESCRIPTION: VPK16 fixes and packs four floating-point numbers from 
vector A into 16-bit quarter words in a single word of 
vector C, packing an array of positive integers with 
values from 8 to 65535, or an array of signed two's 
complement integers with values from -32768 to 32767, 
but does not check for out-of-range values. 


EXAMPLE: 


t 
tw 


i er 
t 


= g 


A. -823 “729 625: 556. 4.1: 324 6255: Led 


Q 


GOB8SFFF8GGG6GSG6 BOG40OO3FFFEDIO1 


A: 8.3 -7.9 6.5 5.6 4.1 3.4 -2.5 1.1 


C : GOB8FFFIGGGEGGGS BUG4VGO3FFFESGIG1 
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* VPKI32 * -——- VECTOR 32-BIT INTEGER PACK --—— * YVPKI32 * 
2 2 * * 
REEKEREKER REEARERREEKE 
PURPOSE: To pack each two 32 bit halfword integer source words 
- into one destination word as halfword-integers-packed. 


CALL FORMAT: CALL VPKI32(A,1I,C,K,N) 


PARAMETERS : A = Halfword integer input vector 
I = Integer element step for A 
C = Halfword~-integer-packed output vector 
K = Integer element step for C 
N = Integer element count (destination words) 
DESCRIPTION: C(m) bits 8 to 31 = A(2m-1) bits 32 to 63 
C(m) bits 32 to 63 = A(2m) bits 32 to 63 


for m=1 to N 
(Bits are numbered 9-63 from left to right). 


VPKI32 packs two halfword integers from vector A into 
32-bit halfwords in a single word of vector C. It 
packs an array of positive integers with values from 
@ to 4294967295, or an array of signed 2's complement 
integers with values from ~-2147483648 to 2147483647. 
VPKI32 does not check for values out of range. 


EXAMPLE: N = 3 
I = 2 
K = 3 (XXX indicates 'undefined') 


A: 8gceeseeadaeaes6 C: KHesssageEgATeoaGA 


8SPABZIESASAIGIASS5 XXXXXXXXXXXXXXXX 
SAEASAGASAGIAST4 XXXXXXXXXXXXXXXX 
8HEGTIBSGISIAGH3 GOBPSIID2GASIS AID 
GISSISSSSSISGSSG2 XXXXXXXXXXXXXXXX 
BIGSISBSASSHSGA1 XXXXXXXXXXXXXXXX 
SLASGIGASAGIOAIGD FFFFFFFEFFFFFFFC 
7FFFFFFFFFFFFFFE 

7FDFFFFFFFFFFFFE 

GOOSIIOSFFFFFFFD 

GOGSSGOSFEEFFFFFC 
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* VSCALE * 
2 * 
RREREEEAEE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


APPENDIX A 


RRKKXRERERE 
x 2 
-—— VECTOR SCALE AND FIX -— * VSCALE * 
x ® 
REEARRRARRER 


To scale the elements of a vector by a power of 2 such 
that a selected scalar will just fit into a specified 


integer bit width, and then fix the scaled elements 
to integers. 


CALL VSCALE(A,1I,B,C,K,N,NB, IEXP) 


= Floating-point input vector 
= Element step for A 
= Floating-point input scalar 


= Element step for C 
= Element count 
= Long-integer input scalar 
(Desired width, 2 to 28 bits, of integers) 
IEXP = Long-integer output scalar . 
(Exponent of scale factor used) 


ZARA WH YP 
a \ 


C(m) = FIX (A(m)*(2**IEXP}) for m=@ to N-1 

where IEXP=NB-E-1, 

and B = FRAC*(2**E). 

VSCALE scales by a power of 2 every element of the 
vector A so that the scalar B will just fit into an 
NB-bit width integer, and then fixes the scaled elements 
and stores them in vector C. IEXP is set to the scale 
factor chosen. If the scalar is larger in magnitude 


than any element of A, no fixing overflows will occur. 


EXAMPLE: 


(with N=5, NB=12) 


B: 18.9 
A: 19.8 5.8 @.2 -4.8 9.91 
C : 1289 649 29° =512 1 
IEXP : 7 
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RRBREEREEERE RERREEKERE 
& & & & 
* VSHFX * ——— VECTOR SHIFT AND FIX -— * VSHFX * 
2 * 4 & 
RERREEEREE RREAEKEERAEE 
PURPOSE: To shift (multiply by a power of 2) and then 


fix (truncate) to integers the elements of a 


= aad ss — Meee te Set Se 


floating-point vector. 
CALL FORMAT: CALL VSHFX(A,1I,C,K,N,NS) 
PARAMETERS : A = Floating-point input vector 
I = Integer element step for A 
C = Long-integer output vector 
K = Integer element step for C 
N = Integer element count 
NS = Integer power of 2 (May be negative) 


DESCRIPTION: C(m)=FIX{A(m)*(2**NS)}; for m=l to N 


EXAMPLE: 
N = 3 
NS = 2 
At 1.0 2.8 3.2 
c 2.4 8 12 
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REEARAKETKEEE RAREKEEKEEE 
® & * *« 
* \vUP8 * -—~- VECTOR 8-BIT BYTE UNPACK --- * VUP8 * 
* c J . = x 
RERELREREE REARREKKKAEE 
PURPOSE: To unpack eight 8-bit unsigned bytes from each 


x - . smh A ms atrinn 
source word and store them in eight destination 


words as 64-bit floating-point numbers. 
CALL FORMAT: CALL VUP8(A,1I,C,K,N) 
PARAMETERS : = Unsigned-byte~integer input vector 
= Integer element step for A 
Floating-point output vector 


= Integer element step for C 
= Integer element count (source words) 


ZROH pY 
u 


DESCRIPTION: Unpacks eight 8-bit bytes from a single word of 
vector A storing them as eight floating-point numbers 
in vector C. The unpacked bytes have values from g to 
253-6 


EXAMPLE: 


Zz 
i] 
N 


Hi 
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REKEEREEEE RREKERKEEK 
& & . 2 . 

* VYUP32 * -—-— VECTOR 32-BIT BYTE UNPACK --—— * VYUP32 * 

> ®& x *& 
RREERKEEREE BEKREREEELR 
PURPOSE: To unpack two 32-bit unsigned halfwords 


° 
from each source word and store them in 


two destination words as 64-bit floating-point 
positive numbers. 


CALL FORMAT: CALL VUP32(A,1I,C,K,N) 


PARAMETERS : Unsigned-halfword-integer input vector 
Integer element step for A 
Floating-point output vector 

Integer element step for C 

Integer element count (source words) 


2A AH DB 
rT) 


DESCRIPTION: VUP32 unpacks two 32-bit halfwords from a single 
word of vector A, storing them as two positive 64-bit 
floating-point integers in vector C. The unpacked 
halfwords have values from @ to 4294967295. 


EXAMPLE: 
N= 4 
A : SIBAISG8SAPSTIAI7 
GISSSSHESDSIGIIS 
GIBSSIS4GSISSSIG3 
GIBSSSS2GAGAISGA1 


C:8.8 7.8 6.8 5.8 4.8 3.9 2.8 1.8 
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* VUPS8 * -——- VECTOR 8-BIT SIGNED BYTE UNPACK ——— * vUPS8 * 
& x & & 
REREKKEEAE PEELE SLE LE £ | 
PURPOSE: To unpack eight 8-bit signed bytes from 


a ee a : _ 


@€ach source-word and store them in eight destination 
words as 64-bit floating-point numbers. 

CALL FORMAT: CALL VUPS8(A,I1I,C,K,N) 

PARAMETERS : A = Signed-byte-integer input vector 

I = Integer element step for A 

C = Floating-point output vector 

K = Integer element step for C 

N = Integer element count (source words) 


DESCRIPTION: VUPS8 unpacks eight 8-bit signed bytes from 
a single word of vector A, storing them as eight 
floating-point numbers in vector C. The unpacked 
bytes have values from -128 to 127. 


EXAMPLE: 


2 
iT] 
N 
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RAREREERERE RREEKEKEEKEE 
= * & & 
* VUPS32 * -—-——- VECTOR 32-BIT SIGNED BYTE UNPACK -—— * YUPS32 * 
* * x * 
RREKREERERE REEARKKEERKEE 
PURPOSE: To unpack two 32-bit signed two's complement 


halfwords from each source word and store them 
in two destination words as signed 64-bit 
floating-point numbers. 

CALL FORMAT: CALL VUPS32(A,1I,C,K,N) 

PARAMETERS : = Signed-halfword-integer input vector 

= Integer element step for A 

Floating-point output vector 

= Integer element step for C 

= Integer element count (source words) 


ZRAQH YP 
rf) 


DESCRIPTION: VUPS32 unpacks two 32-bit signed two's complement 
halfwords from a single word of vector A, storing 
them as two floating-point numbers in vector C. The 
unpacked halfwords have values from -2147483648 to 
2147483647. 


EXAMPLE: 


zm 
iN 


4 

A : GOOGOOS8FFFFFFFI 
SBIsGsseggadggss 
PFFFFFFCFFFFFFFD 
FFFFFFFESS99 9991 


C: 8.9 -7.8 6.09 5.9 =-4.9 <-3.0 -2.2 1.9 
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2 & & x 
* YUUL32 * ——— VECTOR 32-BIT UNSIGNED UNPACK -——— * YUUI32 * 
2 = * x 
REEREREEETE RREKARKEREL 
PURPOSE: To unpack two 32-bit halfword integers from each 


CALL FORMAT: 


source word and store them as two destination words, 
in unsigned integer format. 


CALL VUUI32(A,I,C,K,N) 
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PARAMETERS: A = Halfword integer packed input vector 
I = A address increment 
C = 32 bit integer output vector 
K = C address increment 
N = Integer element count (source words) 
DESCRIPTION: C(2m-1) = A(m) bits g to 31 
C( 2m) = A(m) bits 32 to 63 
for m=9 to N-l 
(Bits are numbered 9-63 from left to right). 
VUUI32 unpacks two 32-bit unsigned halfword integers 
from a single word of vector A and stores them as 
two unsigned halfword integers in vector C. The 
unpacked halfwords have values from 9 to 4294967295. 
EXAMPLE: N = 3 
I = 3 
K = 2 (XXX indicates ‘undefined’ ) 
A: 9@gg9gcsgggegGG7 C: Kesssgsgesgsggdss 
SISISSPGEAAGAAAASS XXXXXXXXXXXXXXXX 
GIGIBAS 4GAGTIGG3 SIDISISIAGAIGI7 
SIVSSSI2GGASFGG1 XXXXXXXXXXXXXXXX 
FRFFFFFFFFPEFEFFE SIgadgagagagggg2 
PFFFFFFDFFFFFFFC XXXXXXXXXXXXXXXX 
FFFFFFFBFFFFFFFA GSIGSSSISSSPSSFS1 
XXXXXXXXXXXXXXXX 
OSGSOSIOFFFFFFFB 
XXXXXXXXXXXXXXXX 
GOOGSGOOFFFFFFFA 
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RREERARETREE RREKEREKER 
4 x ; z & 
* DADD * —-- DOUBLE TO DOUBLE-PRECISION ADD —— * DADD * 
& 2 2 ® 
ERAEREREKEEE RREEKREREE 
PURPOSE: To form a double-precision sum of two 

Anmwhlaeneramiain numhar 

GOUSLE “precision numoers. 


CALL FORMAT: CALL DADD(XDBLE,YDBLE, ZDBLE) 


PARAMETERS : XDBLE = Real vector input (double precision) 
YDBLE = Real vector input (double precision) 
ZDBLE = Real vector output (double precision) 
(A double-precision value is stored ina 
2-element real array. First element contains 
high word, second element contains low word.) 


DESCRIPTION: Adds the double-precision number in XDBLE to the 
double-precision number in YDBLE and stores the high 
word of the double-precision sum in ZDBLE(1) and the 
low word in ZDBLE(2). 
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| PURPOSE: 
| 
| 


{CALL FORMAT: 


DESCRIPTION: 


APPENDIX A 


Perce e ELE ES 
& & 
—-—— DOUBLE ACCUMULATE DOT PRODUCT —-—- * DADOT * 
* * 
RERRERERER 


To perform the dot product of two real vectors, 


4 » s e ° ° ° 
accumulating the result in double precision (128 bits}, 


and returning the result in single precision (64 bits). 


SW = DADOT(N,A,1,B,J) 


Integer element count 

= Real input vector 

Integer element step for A 
= Real input vector 

= Integer element step for 
Real output result 


Lv 


= SUM(A(m) * B(m)) for m=1toN 

= 9.8 for N< 1 
If the element increment, INC, of a vector is negative, 
then the vector is indexed in reverse order, i.e 
element (N-1) * INC + 1 to the first element (BLAS 
convention). 


aD Bh ee 
u 
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* DMUL * 
x & 
RAREKRKERERE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


APPENDIX A 


kkeRKEREEE 
* * 


-——- DOUBLE TO DOUBLE-PRECISION MULTIPLY -—— * DMUL * 


& ® 
RERREREAEE 


To form a double-precision product of two 
double-precision numbers. 


CALL DMUL(XDBLE, YDBLE, ZDBLE) 


XDBLE = Real vector input (double precision) 
YDBLE = Real vector input (double precision) 
ZDBLE = Real vector output (double precision) 


(A double-precision value is stored ina 
2-element real array. First element contains 


nigh w 


| esa2en4 1 - 
ord, secona eiement cont n 


ains low wo 
Multiplies the double-precision number in XDBLE by the 
double-precision number in YDBLE and stores the high 
word of the double-precision product in ZDBLE(1) and 
the low word in ZDBLE(2). 
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REEEERAERER REKREREREE 
2 & & & 
* DNEG * ——-~ NEGATE DOUBLE-PRECISION NUMBER --—— * DNEG * 
; 2 ed & . 
RREEKKEEEX KRERAREKKAKKAE 
PURPOSE: To negate a double-precision number. 


CALL FORMAT: CALL DNEG(XDBLE, ZDBLE) 


PARAMETERS : XDBLE Real vector input (double precision) 
ZDBLE = Real vector output (double precision) 
(A double-precision value is stored ina 
2-element real array. First element contains 
high word, second element contains low word.) 


DESCRIPTION: Negates the doubie-precision number in XDBLE and stores 
the high word of the double-precision result in ZDBLE(1) 
and the low word in ZDBLE(2). 
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* DSUBRR * 
* P 7 
REKKKKKEEE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


APPENDIX A 


RRRERREREE 
* * 


-—-—~— SINGLE TO DOUBLE-PRECISION SUBTRACT -——— * DSUBRR * 


cf & 
REEREKRERE 


To form a double-precision difference of two 
single-precision numbers. 


CALL DSUBRR(X,Y,ZDBLE) 


xX = Real scalar input 
y = Real scalar input 
ZDBLE = Real vector output (double precision) 


(A double-precision value is stored ina 
2-element real array. F 

1igh word, second element contains low word.) 
Subtracts the single-precision number in Y from the 
single-precision number in X and stores the high word of 
the double-precision difference in ZDBLE(1) and the low 
word in ZDBLE(2). 
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RETRETEEEE PEELE SELL £ 
x x & 
* ABS * ——— REAL NUMBER ABSOLUTE VALUE ——- * ABS * 
® 2 & 
RAEETEREKEE RREEEREEEE 
PURPOSE: To compute the absolute value of a real number. 


CALL FORMAT: Function-value 


PARAMETERS: Function-value 
Arg 


DESCRIPTION: Function-value 


FPS 868-7482-881C 


ABS(arg) 


Real Floating-point scalar output 
Real Floating-point scalar input 


larg| 
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RREETEARAEE REEEEEEEER 
& & x x 
* AINT * ——~ TRUNCATE REAL NUMBER -—— * AINT * 
* ® : x x 
KRERREERIEER RRERRAREKE 
PURPOSE: To truncate a real number. 


CALL FORMAT: Function-value = AINT(arg) 


PARAMETERS : Function-value = Real floating-point scalar output 
Arg Real floating-point scalar input 


DESCRIPTION: Function-value = FLOAT(FIXT(arg) ) 
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REKEAREREEE 7 RERRAEEEKRKE 
bi : * 7 & 
* ALOG1G * ——- REAL NUMBER LOGARITHM -—— * ALOGIG * 
. ® & * 
RERERREEELX RRREEEEEKRE 
PURPOSE: To compute the logarithm of a real number. 


CALL FORMAT: Function-value = ALOG(arg) or ALOG19(arg) 


PARAMETERS : Function-value = Real Floating-point scalar output 
Arg Real Floating-point scalar input 


DESCRIPTION: Function-value = Ln(arg); for ALOG 
Log(1%)(arg); for ALOG19 
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REKRERERER REEREREERE 
x & rd & 
* ANINE * -—— ROUND REAL NUMBER TO NEAREST WHOLE -——— * ANINT * 
z & x * 
REREREREEE R&R ReEAARESE 
PURPOSE: To round a real number to the nearest whole number. 


CALL FORMAT: Function-value = ANINT(arg) 


PARAMETERS : Function-value = Real floating-point scalar output 
Arg Real floating-point scalar input 


DESCRIPTION: Function-value = FLOAT(FIX(arg) ) 
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RERRAEEEEKEX REEKEKTAREKE 
bd * & *& 
* ATAN * —— ARCTANGENT OF REAL NUMBER ——— * ATAN * 
& x * * 
RREKREEEEE RRERREEAEEKS 
PURPOSE: To compute the arctangent of a real number 

or of the ratio cf two real numbers. 


CALL FORMAT: Function-value = ATAN(argl) or ATAN2(argl,arg2) 


PARAMETERS : Function-value = Real Floating-point scalar output 
Argl Real Floating-point scalar input 
Arg2 = Real Floating-point scalar input 


DESCRIPTION: Function-value = ATAN(argl) or ATAN(argl/arg2) 
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* CABS * 


REREERERER 


PURPOSE: 


CALL FORMAT: 


PARAMETERS 


DESCRIPTION: 


APPENDIX A 


RERAREREKEE 


——— COMPLEX NUMBER ABSOLUTE VALUE --~- * CABS * 


To compute the 


er te ed 
LLUMVS L eo 


Function-value 


Function-value 
Arg 


Function-value 
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RERREEEKERE 


absolute value (magnitude) of a complex 


= CABS(arg) 


= Floating-point scalar output 
Complex floating scalar input 


SORT (R(arg)**2+I (arg) **2) 
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KRERERERENK RRERRERAELR 
. J * x 
* cCDIV * -———- COMPLEX/COMPLEX DIVIDE --- * CDIV * 
® & & * 
RARRERERERE RRERREEAKEE 
PURPOSE: To divide a complex number into a complex number. 


CALL FORMAT: Function Value = Arg2/Argl 


PARAMETERS : Function Value = Complex Floating scalar output 
Argl Complex Floating scalar input 
Arg2 Complex Floating scalar input 


{R(arg2)+1I(arg2)}/{R(argl)+I(argl) } 


DESCRIPTION: Function Value 
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REREEEEEEE RREREREEEE 
x 2 ® p 
* CDIVRC * —-—- COMPLEX/REAL DIVIDE -——— * CDIVRC * 
, b x ® 
RREKAREREEE REEERAERERE 
PURPOSE: To divide a real number into a complex number. 


CALL FORMAT: Function Value = Arg2/Argl 


PARAMETER: Function Value = Complex Floating scalar output 
Argl = Real Floating-point scalar input 
Arg2 = Complex Floating scalar input 


DESCRIPTION: Function Value = R(arg(2))+I(arg(2))/arg(1l) 
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& & : x x 
* CLOG #* -~—~— COMPLEX NUMBER LOGARITHM --—— * CLOG * 
x & 4 & 
REEEEAERER RERRKERAEREE 
PURPOSE: To compute the natural logarithm of a complex number. 
CALL FORMAT: Function-value = CLOG(arg) 
PARAMETERS : Funct ion-value = Complex floating scalar output 

Arg = Complex floating scalar input 


DESCRIPTION: R(Function-value) ALOG ( (CABS (arg) ) 
I(Function-value) = ATAN(I(arg)/R(arg) ) 
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RAREARREEKEN REREKREEEEEK 
= 5 2 . J 
* CONIG * ——— CONJUGATE OF COMPLEX NUMBER —— * CONJG * 
& ® bd * 
RERREREEEE RELRREREEREE 
PURPOSE: To compute the conjugate of a complex number. 


CALL FORMAT: Function-value = CONJG(arg) 


PARAMETERS : Function-value = Complex floating scalar output 
Arg Complex floating scalar input 


DESCRIPTION: Function-value = R(arg)-I(arg) 
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RERREREEEE REEREREEEE 
5 * & x 
* COSH * ~—— REAL NUMBER HYPERBOLIC COSINE -— * COSH * 
ed & x *& 
REEKETREKEEE RRERETEKEEE 
PURPOSE: To compute the hyperbolic sine or cosine of a real 


oe ee le ee oe 
TIMMS L oe 


CALL FORMAT: Function-value = SINH(arg) or COSH(arg) 


Function-value = Real Floating-point scalar output 
Arg Real Floating-point scalar input 


PARAMETERS : 


DESCRIPTION: Function-value = SINH(arg) or COSH(arg) 
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* CPOWCI * 
* & 
REEKEEAAREE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


--- COMPLEX TO INTEGER POWER -—— 


To raise 
Function 
Function 
Argl 
Arg2 


Function 
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RRAEAKEEEEE 


& 


* 


* CPOWCI * 


& 


x 


REREKREEEEE 


a complex number to an integer power. 


Value 


Value 


Value 


Argl**Arg2 


Complex Floating scalar output 
Complex Floating scalar input 


Integer scalar input 


{R(argl)+I(argl) }**arg2 
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* 2 : & . 
* CPOWRC * ——— REAL TO COMPLEX POWER -——— * CPOWRC * 
® & ® * 
RERERARAEEE REREREREEE 
PURPOSE: To raise a real number to a complex power. 


CALL FORMAT Function Value = Argl**Arg2 


PARAMETERS : Function Value = Complex Floating scalar output 
Argl Real Floating-point scalar input 
Arg2 = Complex Floating scalar input 


DESCRIPTION: Function Value = argl**(R(arg2)+I(arg2) ) 
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RERREKEURR RREREREEKE 
 d 7 x 
* csoRT * —— SQUARE ROOT OF COMPLEX NUMBER -——— * CSQRT * 
*& b x * 
REREKEEREE REREAREREEE 
PURPOSE: To compute the square root of a complex number. 


CALL FORMAT: Function-value = CSQRT(arg) 


PARAMETERS : Function-value = Complex floating scalar output 
Arg Complex floating scalar input 


DESCRIPTION: if R(arg) > @ R( function value) = F 
I(function value) = I(arg)/(2*F) 
if R(arg) < g R(function value) = I(arg)/(2*F) 
I(function vaiue) = SIGN(I(arg))*F 


where F = SQRT((ABS(R(arg) )+CABS(arg))/2) 
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= EXP * 


RAERERAAAE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


APPENDTX A 
kekERREREKE 
* te 
——- EXPONENTIAL OF REAL NUMBER —— * EXP +* 
* # 
keeRREREEEX 


To compute the exponential of a real number. 


Function-value = EXP(arg) 


Function-value = Real Floating-point scalar output 
Arg Real Floating-point scalar input 


Function-value = Exp(arg) 


NOTE: arg>799.889 traps with an overflow error 
condition. 
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REEEREREREE xz 
& *& x 
* IDIM * ——- INTEGER/INTEGER POSITIVE DIFFERENCE ——— * 
* Pd . J 
RRRAEREREEE zak 
PURPOSE: | To compute the integer positive difference of two 


integers. 
CALL FORMAT: Function-value 
PARAMETERS : Function-value 
Argl 
Arg2 


DESCRIPTION: Function-value 


= IDIM(argl,arg2) 
= Integer scalar output 
= Integer scalar input 


= Integer scalar input 


= MAX( (argl-arg2),@) 


RREREERE 
& 
IDIM * 
& 
SERUREEE 
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REKKAERERE RREREREREE 
& _® *& 
* Ipow * —— INTEGER TO INTEGER POWER -—- * Ipow * 
* & * & 
RRREKEEEEE REKKREKEREE 
PURPOSE: To raise an integer number to an integer power. 


CALL FORMAT: Function Value = Argl**Arg2 
PARAMETERS : Function Value = Integer scalar output 
Argl = Integer scalar input 


Arg2 = Integer scalar input 


DESCRIPTION: Function Value = argl**arg2 
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* MOD +* 
® x 
RERRREEEEE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


——~— INTEGER/INTEGER DIVIDE REMAINDER --—— 


To compute the 
by another. 


Function-value 
Function-value 
Argl 
Arg2 


Function-value 
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remainder when one integer 


MOD(argl,arg2) 
Integer scalar output 
Integer scalar input 


Integer scalar input 


Argl-INT(argl/arg2)*arg2 


RRREREREEK 
& & 


* MOD * 
* * 
RERRKEKRERE 


is divided 
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* RAN * 


RRREEREREE 


PURPOSE: 


CALL FORMAT: 


PARAMETERS : 


DESCRIPTION: 


EXAMPLE: 
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RREKRERETEE 
4 £ 
—-- SCALAR RANDOM NUMBER GENERATOR ——- * RAN # 
* &® 
PREECE LE ST 


To generate one pseudo-random number. 
Function-value = RAN(SEED) 


Function-value = Floating-point output scalar 
Output random number 

SEED = Integer input/output scalar 
Input: random number seed 
Output: last integer generated 


RAN returns one pseudo-random fioating-point 
number between 9.8 and 1.9. The routine uses 

a linear congruential generator initialized by 
SEED to generate an integer, which is then scaled 
to produce the function-value. SEED is replaced 
with the integer generated. SEED may be any 
integer between 9 and 2**26-1. 


RAN generates the same sequence of integers as VRAND. 
Thus the two statements 

C = RAN(SEED) 
and 


co 


are equivalen 


SEED = 1999 


RAN(SEED): @.8994849494996693 
SEED : 53719635 


Ey 
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KRREREEREEE . RRRRAREKEER 
x & . & 
* RPOW * ——— REAL TO REAL POWER ——— * RPOW * 
x *® & x 
RREKERERAER REEKEREKEER 
PURPOSE: To raise a non-negative real number to a real power. 


CALL FORMAT: Function Value = Argl**Arg2 


PARAMETERS : Function Value = Real Floating-point scalar output 
Argl Real Floating-point scalar input 
Arg2 Real Floating-point scalar input 


DESCRIPTION: Function Value = argl**arg2 


whole number, Argi can be negative.) 


wi 
fu 


(if ArgZ i 
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RERRRAEAEUE RREEREEREE 

b & & & 

* RRCP * ——- REAL RECIPROCAL -—- * RRCP # 

® & & x 
REEREREEEE REEKKKERKKE - 
PURPOSE: To divide a real number into a real number or into l. 


CALL FORMAT: Function Value = Arg2/Argl 
or 1.9/Argl 


PARAMETERS : Function Value = Real Floating-point scalar output 
Argl = Real Floating-point scalar input 
Arg2 = Real Floating-point scalar input 


DESCRIPTION: Function Value = arg2/argl for RDIV 


ae 1 Os 1 #¢£ 
Or 4.6/a4rG.i itor RRCP 
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RRETREREEX RRKERKEEEE 
* * * * 
* SIGN * ——- REAL NUMBER SIGN TRANSFER -——- * SIGN * 
® 3 & & 
RREERARREREE RERREREKEE 
PURPOSE: To give the magnitude of a real number with the sign 
of a second real number. 
CALL FORMAT: Function-value = SIGN(argl,arg2) 
PARAMETERS : Function-value = Real Floating-point scalar output 
Argl = Real Floating-point scalar input 
Arg2 _ = Real Floating-point scalar input 
DESCRIPTION: Function-value = Sign{arg2)*ABS(argl) 
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* & 


* SINCOS * 
* * 


RRAEREKKERE 


PURPOSE: 


CALL FORMAT: 


REREKERARER 
& *& 


-——- REAL SINE AND COSINE --—— * SINCOS * 
* *. 


RAEREKERERE 


To compute the sine and cosine of a real number. 


CALL SINCOS(A,CA,SA) 


PARAMETERS : A = Floating-point input scalar 

CA = Floating-point output scalar 

SA = Floating-point output scalar 

DESCRIPTION: CA = COS(A) 

SA = SIN(A) 

SINCOS computes both the sine and the cosine in 

about the same time as the SIN function alone. 

NOTE: A 32-bit integer overflow exception is generated 
if the input argument is too large (greater 
than approximately 8.9E+5). In this case, the 
output result has less than six decimal digits of 
precision. 

An added feature of this routine is that it can also 

be called as a complex function. If FIFSPR_SINCOS 

is declared as complex, the call 

Function-value = FIFSPR_SINCOS(A) 
returns the complex value 
Function-value = CMPLX(COS(A),SIN(A)). 
This is convenient for converting polar coordinates 
to rectangular coordinates. 
EXAMPLE: 


A = 9.8 
CA = 1.9 
SA = 8.2 
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RETEREREELX 
& ® 


* SORT * 
* * 


RREREEEEKEE 


PURPOSE: 
CALL FORMAT: 


PARAMETERS: 


DESCRIPTION: 
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RERKRAEKEKER 
& * 


-—— SQUARE ROOT OF REAL NUMBER --—— * SQRT * 


To compute the 
Function-value 


Function-value 
Arg 


Function-value 


& x 
REEKEEEERX 


Square root of a real number. 


SORT(arg) 


= Real Floating-point scalar output 
Real Floating-point scalar input 


= SQRT(arg) 
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RREARRREREETE SRRRREKTEEX 
e * ey & * 
* ‘TANH * --- REAL NUMBER HYPERBOLIC TANGENT -—-—-— * TANH * 
e & & * 
RRREEKAEEAE REEKRERAREEE 
PURPOSE: To compute the hyperbolic tangent of a real number. 


CALL FORMAT: Function-value = TANH(arg) 


PARAMETERS : Function-value = Real Floating-point scalar output 
Arg Real Floating-point scalar input 


DESCRIPTION: Function-value = TANH(arg) 
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‘DATA REPRESENTATIONS FOR STORING 
SPARSE VECTORS AND MATRICES 


B.1 INTRODUCTION 


This appendix presents information to help the user understand and use 
the sparse vector and sparse matrix subroutines. It describes the data 
representations (or formats) both accepted as input and produced as 
output by these routines. This appendix also spells out parameter 
Naming conventions common to many of these subroutines. 


There are four subroutines that convert sparse vectors and matrices 
between their packed and full representations: Sparse Vector Pack 
(SVPACK), Sparse Vector Unpack (SVUPCK), Sparse Matrix Pack (SMPACK), 


and Sparse Matrix Unpack (SMUPCK). 


B.2 SPARSE VECTOR STORAGE 


An N-dimensional sparse vector V is represented in packed-vector format 
by N, NS, S, and IEN where: 


N a scalar, is the dimension of V. 

NS a scalar, is the number of nonzero values in V. 

s a vector of length NS, contains the nonzero values of V. 
IEN a vector of length NS, contains the location in V of each 


corresponding element in § [{i.e., V(IEN(k)) = S(k) for k=1,NS]. 
For example, the following sparse vector 
(8.8 3.2 8.86 7.8 9.8 8.89 8.8 -19.3] 


can be represented in packed-vector format as follows: 


N: 8 
NS: 3 

S: [3.2 7.8 19.3] 
IEN: {2 4 8 j 


So, S(1)‘'s location in V can be found in IEN(1), S(2)'s in IEN(2), ..., 
S(NS)'s in IEN(NS). 


The nonzero values in S are generally ordered as they appear in V. 


However, they can be ordered differently if the order is compatible 
with the subroutine to be used. 
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Except for differences in the IP vector, formats I and III are the 
same, as are formats II and IV. 


Each attribute associated with a particular format type and the 


consequences of using that attribute are explained in detail 
in the sections that follow. 


el Matrix Format Tvne I (COL—ORDER PTRS_ONT.Y) 


A sparse matrix A is represented by M, N, NS, S, IN, and IP, 
in format I where: 


M a scalar, is the number of rows in A. 

N a scalar, is the number of columns in A. 

NS a scalar, is the number of nonzero values in A. 

s a real vector of length NS, contains the nonzero values of A 
in column order. 

IN an integer vector of length NS, contains the row in A of 
each corresponding value in § [i.e., IN(k) = row in A of S(k) 
for k=1,NS]. 

IP an integer vector of length N+tl, contains one element for every 
column in A. 
Each element indicates the location in S that holds that column's 
first nonzero value (exception: empty columns). 
IP's N+tlst element is a sentinel. 


The sentinel element IP(N+1) holds the number NS+1. 
In general, IP(i) contains the location in S that refers to A's 


Sek h FPATiIuMmn far i 
27CGR COLUR, LVL 4 


If a column in A is empty, then the entry in IP for that column 
is the same as the entry for the next nonempty column, or if 
there is no such column, 

sentinel value in IP(N+l) is used. 

The matrix: 


9.9 
0.8 
9.2 
G.8 
9.8 


as expressed in Type I Format: 


Ms: 5 
N: 6 
NS: 8 
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B.3.3 Matrix Format Type III (COL-ORDER PTRS—SUMS) 


A sparse matrix A is represented by M, N, NS, S, IN, and IP, in format 
III where: 


M a scalar, is the number of rows in A. 

N a scalar, is the number of columns in A. 

NS a scalar, is the number of nonzero values in A. 

S a real vector of length NS, contains the nonzero values of A in 


column order. 


IN an integer vector of length NS, contains the row in A of each 
corresponding value in S [{i.e., IN(k) = row in A of S(k) for 
k=1,NS]. 

IP an integer vector of length 2*N, contains two elements for every 
column in As: 

(a) the location in S that holds that column's first 
nonzero value (exception: an empty column). 
(b) that column's total number of nonzero elements. 


IP(i) and IP(i+N) always refer to the i-th column in A, for i=1,N. 
IP(1) to IP(N) holds locations as in (a) above and IP(N+l) to IP(2*N) 
holds sums as in (b) above. 


If a column in A is empty, then the (a)-entry in IP for that column is 
the same as the (a)-entry for the next nonempty column, or if there is 
no such column, the number NSt+l. (Note that the (b)-entry is zero.) 


The matrix: 9.9 G€.8 6.89 G€.89 O.8 G.8 
$9.08 4.5 G.8 9.2 3.8 9.f 
9.9 8.89 G€.89 €.8 G.8 GF 
8.8 9.9 7.1 §.8 3.8 G.8 
9.8 1.3 G.82 8.3 G.8 G.G 
as expressed in Type III Format: 
M: 5 
N: 6 
NS: 8 
S: 14.5 9.9 1.3 7.1 9.2 5.8 8.3 3.8] 
IN: [2 4 5 4 2 4 5 2 ] 
IP: £1 1 4 5 8 9 g 3 1 3 1 9) 


Note that lengths of S and IN equal NS (=8); S is in column order; the 
length of IP equals 2*N (=12); IP contains both locations and sums; IN 
contains row numbers. 
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SPARSE LINEAR SYSTEM ROUTINES 


C.l1 INTRODUCTION 
Thi ndix contains information ¢: 
the sparse linear system routines in the Advanced Math Pibesey: ‘The 
Sparse linear system routines are APAL64 routines that provide an 
efficient method for solving the linear system Ax = b where the 
coefficient matrix is sparse and is stored in packed form. 


There are twelve generic sparse linear system routines in all. The 
name of each routine consists of a four-letter generic name followed by 
the single digit "2". The first two letters of the name indicate the 
coefficient matrix type (i.e., the problem domain), and the last two 
letters indicate its function. The single digit is a version number 
and is not included on the names of the original routines, which were 
superseded as of the F93 release (see Appendix G). 


The types of coefficient matrices are: 


RU A is real. 

RS A is real and symmetric. 
CU A is complex. 

cs A is complex and symmetric. 


The functions performed are: 
FR Factor the coefficient matrix. 


SV Solve the system given the factorization of the 
coefficient matrix. 


FS ’ Factor and solve (combines FR and SV). 


In general, the time required to factor the coefficient matrix is much 
greater than the time required to solve the factored system. 

Therefore, by having separate routines for each of these functions, the 
factorization need only be performed once when solving a number of 
systems that all have the same coefficient matrix. 
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Denote the determinant of a square matrix A by Det(A). The “not equal" 
relation will be denoted by the symbol "#". 


Assume an n x n lower-triangular matrix L, and n x n upper-triangular 
matrix U, such that A = LU. Then the system Ax = b is equivalent to 
LUx = b. Letting Ux = y, where y is an n-dimensional vector, then the 
system becomes Ly = b. Thus, it is possible to decompose the original 
system into two triangular systems which, in general, are easier to 
solve. It is then possible to find the solution to the original system 
x, by the following two steps: 


1) Solve Ly = b for y by forward elimination 


2) Solve Ux = y for x by backward substitution 


_ If there does exist an L and U such that LU = A, then L and U are not 
uniquely determined unless additional conditions are imposed. One such 
set of conditions is to require the following: 


U(i,i) = 1 for i=lton. 


By imposing this restriction on U, the remaining elements of L and U 
can now be solved obtaining the following: 


A(i,j) - Sum(L(i,k) * U(k,j), k=1,4-1]) 
for i=1ton, j = 1 ton, and i >= j eq(la) 


GB(i,j) 


U(i,j) = (ACi,j) - Sum(L(i,k) * U(k,j), k=1,i-1]) / L(i,i) 
for i=1ton-l, j = 2 ton, and j >i eq(1b) 


It is clear from an examination of the expressions above that a unique 
L and U exist if and only if L(i,i) # @ for i=1ton-l. Letting A{k} 
denote the k-th order principle submatrix of A (i.e., the submatrix 
formed by the intersection of the first k rows and the first k columns 
of A), then it follows from equation (1) that A{k} = L{k}u{k}. Recall 
from elementary linear algebra that: 


(a) if A = BC, then Det(A) = Det(B)Det(C); and 
(b) if T is ann x n triangular matrix, 


then Det(T) = Prod({T(i,i}, i=l,n]}. 
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A common variation of the method of LU factorization involves the 
further factorization of L into MD where M is a lower-triangular matrix 
with M(i,i) = 1 for i=l to n and D is a diagonal matrix. The elements 
of M and D are found to be: 


M(i,j) = L(i,3) / L(j,j) eq(2a) 
D(i,i) = L(i,i) eq(2b) 


tions (1) and (2 


a d can be used to show that M is the transpose of J 
A is symmetric. 


\ 
i 
The LDU theorem can now be stated. 


C.3.2 LDU Theorem 


If A is ann x n matrix, then there exist unique matrices L, D, and U, 
where L is lower-triangular with L(i,i) = 1, D is diagonal with D(i,i) 
# 9, and U is upper-triangular with U(i,i) = 1 such that A = LDU if and 
only if Det(A{k}) # @ for k = 1 to n. Furthermore, if A = LDU and A is 
symmetric, then L is the transpose of U. 


If A is factored into LDU, then the original system, Ax = b, is 
equivalent to LDUx = b. Letting Ux = y and Dy = z where y and z are 
n-dimensional vectors, then the original system decomposes into two 
triangular systems and a diagonal system that are solved by the 
following three steps: 


1) Solve Lz = b for z by forward elimination. 
2) Solve Dy = z for y. 
3) Solve Ux = y for x by backward substitution. 


Since LDU-factorization requires more work than LU-factorization, the 
later is preferable unless A is symmetric. In that case, the direct 
computation and storage of U is unnecessary since U is the transpose of 
L and the factors are written LDL’. 


C.4 FILL-IN 


If the coefficient matrix A is sparse, (this is assumed when using the 
sparse system routines) store only the nonzero elements of A with 
information about the location of the nonzero elements. (The manner in 
which this is done is described in Section C.5.) It is very desirable 
to do this since both storage requirements and execution time can be 
greatly reduced. 
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The following algorithm is given in the form of a FORTRAN subroutine 
for determining fill-in: 


SUBROUTINE FILLIN(N, A, IA) 


Cc 
Cc GIVEN AN N BY N MATRIX, A, THIS ROUTINE RETURNS AN N BY 
Cc N LOGICAL MATRIX, IA, WHERE IA(I,J) IS TRUE IF A(I,J) IS 
Cc A SPARSE ELEMENT AND FALSE OTHERWISE 
Cc 

REAL A(N,N) 

LOGICAL IA(N,N) 
Cc 


pO llg@ I=1,N 
TA(I,1) = .FALSE. 
IA(1,I) = .FALSE. 
IF(A(I,1) .NE. 9.8) IA(I,1) 
IF(A(1,I) .NE. 9%.) IA(1,T) 
118 CONTINUE 


- TRUE. 
- TRUE. 


Cc 
DO 158 J = 2, N 
DO 149 I = 2, N 
IF(A(I,J) .NE. 9%.8) GO TO 139 
K2 = MIN@(I,J) -1 
DO 128 K = l, K2 
IF(IA(I,K) .AND. IA(K,J)) GO TO 139 
128 CONTINUE 
IA(I,J) = .FALSE. 
GO TO 149 
138 CONTINUE 
 IA(I,J) = .TRUE. 
14g CONTINUE 
158 CONTINUE 
RETURN 
END 


The amount of fill-in varies as the rows and columns a A are permuted 
and algorithms exist to minimize the fill-in. However, any permuting 
of the rows and columns of A to decrease fill-in may be detrimental to 
the numerical stability. 


Before leaving the subject of fill-in, note that if A is a band matrix, 
then the superposition of L and U will also be a band matrix and will 
have the same bandwidth as A. Therefore, if A is a band matrix where 
the nonzero elements are dense within the band consider every element 
within the band to be sparse without introducing a great number of 
unnecessary sparse elements. 
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Finally, if A if unsymmetric, then an additional integer vector IDP of 
length N is required for pointers into S to the diagonal elements of A. 
For example, 

e If A is real and A(j,j) is stored in S(k), then IDP(j) = k. 


e If A is complex and A(j,j) is stored in S(2*k-1) and S(2*k), 
then IDP(j) = k. 


Consider the following example; let A be the real matrix. 


Note that A(4,5) is a sparse element since it is a fill-in element. 
The vectors S, IRN, ICP, and IDP that are required to represent A are: 


WORD Ss IRN LCP IDP 


o 


UQNH PH We dN 
NMRMQ NN M MM HY 


' 
. ¢ 
Od fF NM 


OMNI AWN WN KH 
Om N Pe WwW NF 


The output from the factorization routines and the input to the 
solution routines require these same vectors except that S then 
contains the sparse elements of the superposition of L and U on A (L', 
D, and U if A is symmetric) with the diagonal elements replaced by 
their reciprocals. (See the example above.) . 


Therefore, the superposition of L and U with the diagonal elements 
replaced with their reciprocals is 
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BASIC LINEAR ALGEBRA SUBPROGRAMS 


D.1 INTRODUCTION 


; mation to help the user understand and use 
the routines, which conned bate the basic linear algebra subprograms 
(BLAS) as implemented within the LINPACK Users' Guide Manual, Appendix 
A. These routines are a subset of the basic linear algebra subprograms 
developed by Lawson, Hanson, Kincaid, and Krogh (refer to ACM Trans. 
Math. Software 5, 3 (Sept. 1979) pp. 324-325) for many of the basic 
vector operations of numerical linear algebra. The package was 
intended to be called from FORTRAN programs, and was developed to focus 
on performance improvements of the well known set of LINPACK routines 
(refer to the LINPACK Users' Guide, Appendix A). 


In addition, four routines have been added which are extensions to four 
of the BLAS routines (real and complex versions of the dot product and 
scalar times vector plus vector) which provide for repeated invocations 
with only one subroutine call. These are useful in many applications 
including matrix multiply and matrix factoring (refer to examples 
D.4.9, D.4.19, and D.4.11). 


Double precision entry points allow the routines to handle standard 
calls to BLAS double-precision routines. There are no specific 
double-precision routines implemented, since the single precision 


ee a es alee ete De wot Ae Soe eee ee aa 


routines use tne standara 64-bit wilge flcoating—- point MUMDeLS. 


When called from FORTRAN, the BLAS routines perform according to the 
algorithmic description in Appendix A, LINPACK User's Guide. In 
particular, negative subscript increment specification results in 
adjustment of the vector base address, as described in Section D.2. 
(No such base address adjustment needs to take place when the MLSP 
entries are used. However, when calling the routines from APAL64 base 
address adjustment is used.) 


Much of the information in Sections D.2 and D.4 is taken from Appendix 
3 of the NTIS-distributed Sandia National Labs. report, SAND77-0898, 
Basic Linear Algebra Subprograms for Fortran Usage, by Lawson, Hanson, 
Kincaid, and Krogh, and is reprinted with their kind permission. 
Floating Point Systems, Inc., gratefully acknowledges the suggestions 
given by R. J. Hanson. 
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D.3 ROUTINE CALLING SEQUENCES, ALGORITHMS, TIMINGS 


The names of entities used in BLAS calls conform in general to standard 
FORTRAN conventions. In particular, names that begin with I or N 
pertain to integer data types; names that begin with C pertain to 
complex data types, and names that begin with S (for scalar) pertain to 
real (floating-point) data types. 


The roots of the names pertain to function. The routines with -DOT- as 
root Calculate different versions of the dot product, SDOT calculating 
the inner product of real vectors, CDOTC and CDOTU calculating complex 
inner products conjugated and unconjugated respectively. 


COPY Replaces (moves or copys) elements of a vector with elements 
of another. 

AXPY Stands for "ax+y". It is intended to perform the elementary 
matrix operation of adding to the elements of a vector the 


scalar multiple of another vector. 


SCAL Multiplies a vector by a scalar. 


SWAP Interchanges (or swaps) elements of two vectors. 


ASUM Calculates the absolute sum of a vector; that is, the sum of 
the absolute values of each element. 


I-AMAX Calculates the index, or subscript, of the component of a 
vector of the largest absolute value. 


S-NRM2 Calculates the 2-norm, or Euclidean length of a vector. It 
carefully concerns itself with scaling problems to maintain 
accuracy and exponent range, by testing each component before 
adding its square to the accumulating partial sum. Usually it 
would be appropriate to use SQRT(DOT) for the same operation 
with greater speed but less robustness. 


ROT Rotates a vector of pairs of points. 


The parameter names are also standardized. These routines all deal 
with one or two vectors, usually coming from matrix rows or columns. 
The first vector is X; the second, -Y. Increments between consecutive 
elements of a vector are named INCX and INCY. Scalars are named 

A and -B. 


Speed values reflect average values, without regard for vector. 
placement, for typical APFTN64 compilations. Often much improvement is 
possible by judicious placement of elements among memory modules. 

Alss, initial setup times are not included, only the loop values, which 
results in a value which is a constant multiple of N, the number of 
elements in the destination vector. 
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D.3.4 Complex Function CDOTU(N,CX, INCX,CY, INCY) 


Function value = sum( CX(m)*CY(m), for the N vector elements 
indexed by m). 


D.3.5 Subroutine CROTG(CA,CB,SC,CSIN) 


sc := |CA|/r, CSIN := conjugate(CB)*CA/|CA|/r, CA := CR 
where: rc=sqrt(|CA|**2 + |CB}**2) and SC,CSIN chosen to satisfy 
CR = SC*CA+CSIN*CB 
g = CSIN'*CA+ SC*CB. 


D.3.6 Subroutine CSCAL(N,CA,CX, INCX) 


CX(m) := CA*CX(m), for the N vector elements indexed by m. 
D.3.7 Subroutine CSSCAL(N,SA,CX, INCX) 
by m. 


CX(m) := SA*CX(m), for the N vector elements indexed 


D.3.8 Subroutine CSROT(N,CX, INCX,CY, INCY,SC,SS) 


CX(m):= SC*CX(m)+SS*CY(m) 
CY(m) :=-SS*CX(m}+SC*CY(m), for the N vector elements indexed by m. 


D.3.9 Subroutine CSWAP(N,CX, INCX,CY, INCY) 


CX(m) :=: CY(m), for the N vector elements indexed by m. 


D.3.18 Integer Function ICAMAX(N,CX, INCX) 


= I such that |Re CX(I)|+|Im CX(I)| is largest of 


Function value = 
the N values |Re Cx(m)|+|Im Cx(m)]|. 


D.3.11 Integer Function ISAMAX(N,SX,INCX) 


= smallest I such that |SX(I)| is largest of all N 


Function value = 
values |SX(m)|. 


D.3.12 Real Function SASUM(N,SX, INCX) 


sum( |SX(m)|, for the N values indexed by m). 


Function value = 
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D.3.21 Subroutine SROTM(N,SX, INCX,SY,INCY, PARAM) 


If PARAM(1) = 1.8 then 


SX(m) 
SY(m) 


= PARAM(2)*SX(m) + SY (m) 
= -SX(m) + PARAM(5)*SY(m), 


for the N vector elements indexed by m. 


If PARAM(1) = 8.8 then 
SX(m) := SX(m) + PARAM(4)*SY(m) 
SY¥(m) := PARAM(3)*SY(m) + SY(m), 


for the N vector elements indexed by m. 
If PARAM(1) = -1.9 then 


SX(m) 
SY(m) 


= PARAM(2)*SX(m) + PARAM(4)*SY(m) 
= PARAM(3)*SY(m) + PARAM(5)*SY(m), 


for the N vector elements indexed by m. 
If PARAM(1) is not 1, 9, or -l, the routine returns without modifying 


the vector elements. It thus becomes equivalent to an identity 
transformation. 


D.3.22 Subroutine SROTMG(D1,D2,Bl1,B2, PARAM) 


PARAM(1) := 9.9 
PARAM(3,4) := -B2/Bl, D2*B2 / D1*Bl, so that the SROTM matrix 
becomes (H11,H21,H12,H22) = (1,-B2/Bl1,D2*B2/D1*Bl,1). 
Dl : 
D2 := D2/U 

:= B1*U where U = 1.9 + (D1*B1*B1)/(D2*B2*B2). 
Tf |D1*B1*B1! =< |pD2*B2*B2!| and D2*B2 <> @ then 
PARAM(1) := 1.98 
PARAM(2,5) := D1*Bl1/(D2*B2) , B1/B2 so that the SROTM matrix 
becomes (H11,H21,H12,H22) = (D1*B1/D2*B2,-1,1,B1/B2). 
D1,D2,Bl := D2/U,D1/U,B2*U where U = 1 + D1*B1*B1/(D2*B2*B2). 
If D2*B2 = 8, then 


the rotation matrix in SROTM becomes the identity, PARAM(1) 
s= -2.9 
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Memory words occupied by X may intersect those occupied by Y. In fact, 
X and Y may coincide. However, memory occupied by Z should not, in 
general, intersect that occupied by X or Y. 


If N < 1, SDOTN returns with no action taken. 


If M < 1 and ISW{1] 1, SDOTN returns with no action taken. 


If M < 1 and Isw[{1] 9, SDOTN returns with Z2(j) = 9.8 for j = 1 toN. 


In general, M < 1 implies a zero sum of products. 


D.3.26 Complex Subroutine CDOTN(ISW,N,M,X,IXI,IX0,Y,1YI,1YO,2,12Z0O) 
Z(jz) = rc * C(jz) + s * SUM[ A(ix) * Bliy), i=1,M] j=1,N 
where: ix = (j-l) * IXO + (i-1) * IxXI +1 

iy = (j-l) * IYO + (i-l) * I¥yI +1 
jz = (j-l) * 120 +1 


s = 1.9, if ISW[9] = g 
= -1.%8, if Isw{[g] = 1 


r= @.8, if ISW[l] = 9 
= 1.8, if Isw(1] = 1 


A = X , if ISW([2] = g 
Conjg(X), if ISW([2] = 1 


B = ¥Y , if ISW[{[3] = ¢g 
= Conjg(Y), if ISW{3] =1 
Cc = Z, if ISW([4] = 9 


Conjg(Z), if ISW([4] = 1 , 


and ISW[(k] bit k of ISW. 

ISW is a one word function selector switch and is treated as a bit 
string with the bits numbered from the least significant bit (bit @). 
If a given bit is set (equal to one), then the function option that 
corresponds to that bit is selected. 


If IZ0 = g, then CDOTN sets 2(1) equal to the accumulated sum of all N 
dot products. If ISW{1] = 1 also, then input Z(1) is added to this 
sum. 


Memory words occupied by X may intersect those occupied by Y. In fact, 
X and Y may coincide. However, memory occupied by Z should not, in 


pe ee ae ee 


i | ~ mou meen 2s a 
enerail, intersect that sccupied 
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D.3.28 Subroutine CAXPYN(ISW,N,M,A,IAO,X,IXI,IxX0O,Y,IYI,1YO) 
Y(iy) = s * B(ja) * Z(ix) + Y(iy), i=1,M j=1,N 
where: ja = (j-1) * IAO +1 

ix = (j-1) * IxO + (i-1) * IxXI +1 
iy = (j-1) * IYO + (i-l1) * IYI +1 


s = 1.9 , if Isw(g] = @ 
= @, if Isw(gj = 1 


B = A, if ISW(2] = 2 
Conjg(A), if ISW([2] = 1 


Z = X , if ISW[3] = g 
Conjg(X), if ISW([3] = 1 


and ISW({k] = bit k of ISW. 


ISW is a one word function selector switch and is treated as a bit 
string with the bits numbered from the least significant bit (bit @). 
If a given bit is set (equal to one), then the function option that 
corresponds to that bit is selected. 


Memory words occupied by A may intersect those occupied by X. However, 
memory occupied by Y should not, in general, intersect that occupied by 
A or X. 


Furthermore, the user will not get meaningful results when distinct 
“columns" of Y¥ intersect. For instance, if M = 19@, IYI = 1 and IYO = 
96, then Y(97,1) = Y¥(1,2), Y(98,1) = Y(2,2) etc. 

However, cases involving IYO = 9 produce meaningful results in that the 
products are accumulated to Y. That is, successive results bound for 
the same storage location in Y are added together rather than stored 
over each other. In this case, the calculation is reduced to a single 
call to CDOTN which executes much faster than the general case speeds 
given in the routine documentation. 

IYI = 8 is of no real value and is omitted for speed and simplicity. 
If N < 1, CAXPYN returns with no action. 

If M < 1, CAXPYN returns with no action. 


If IYI = 9, CAXPYN returns with no action. 
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D.4.5 Set to Identity 


Given an N by N matrix A, to set A = the identity matrix and then 
B= A. 


DO 58 J=1,N 
59 CALL SCOPY(N,%.E9,8,A(1,J),1) 
CALL SCOPY(N,1.E0,0,A,MDA+1) 


DO 69 J=1,N 
rap 4 ao rola tnbh?d Rs 
On CALL SCOPY(N,A{(1,J3),1,B(1,75),1) 


D.4.6 Matrix Columns Interchange 


To interchange the columns of an M by N matrix C, where the column to 
be interchanged with column J is in a type INTEGER array IP(*), and has 
the value IP(J). 


DO 79 J=1,N 

L=IP(J) 

IF(J.NE.L) CALL SSWAP(M,C(1,J),1,C(1,L),1) 
78 CONTINUE 


D.4.7 Matrix Transposition 


To transpose an N by N matrix A in-place, where MDA is the first 
dimensioning parameter of the array A(*,*). 


IF(N.EQ.1) GOTO 85 
DO 69 J=1,N-1 
8g CALL SSWAP(N-J,A(J,J+1),MDA,A(J+1,J),1) 
85 CONTINUE 
D.4.8 Column Vector Circular Shift 


Finally, an inefficient but illustrative code segment which swaps 
in-place the components of the column vector 


(x1,...,XK,xKt1l,...,xXN) 
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D.4.19 Matrix Factorization Using SAXPYN 


This subroutine performs matrix factorization A=LU without pivoting > 
using SAXPYN. L replaces the lower part of A excluding the diagonal. 
L is assumed implicitly to have 1's on its diagonal. U replaces the 
upper part of A including the diagonal. A itself is treated as a 
doubly dimensioned array with first dimension NO. A is assumed to 
contain an NI x NI matrix stored by rows rather than the usual storage 
by columns. This storage scheme allows SAXPYN to more efficiently 


ee es -_ faw a ee 


ee ee ee we sem ~ 
PLrucessa the CULLeUL Liuw bei my “seu vue elimination. 


SUBROUTINE MFBGE(A,NI,NO) 
REAL A(1) 
INTEGER NI,NO 


Cc 
IF(NI.LE.1) RETURN 
JINV=1 
NOP=NO+1 
Cc 
DO 198 I=1,NI-1 
AINV=1.8/A(JINV) 
JC=JINV+NO 
Cc 
Cc COMPUTE THE NEXT COLUMN OF L 
Cc 
CALL VSMUL(A(JC),NO, PEIN A(JC),NO,NI-I) 
MN=NI-I 
Cc 
Cc PERFORM THE ELIMINATION GETTING A NEW LOWER RIGHT MINOR 
Cc 
CALL SAXPYN(1,MN,MN,A(JC),NO,A( JINV+1),1,0,A(JIC+1),1,NO) 
Cc 
JINV=JINV+NOP 
189 CONTINUE 
Cc 
RETURN 
END 


D.4.11 Matrix Factorization Using SDOTN 


This subroutine performs matrix factorization A=LU without pivoting 
using SDOTN. L replaces the lower part of A excluding the diagonal. 
L is assumed implicitly to have 1's on its diagonal. U replaces the 
upper part of A including the diagonal. A itself is treated as a 
doubly dimensioned array with first dimension NO. A is assumed to 
contain an NI x NI matrix stored by columns. Doolittle's method is 
used. 
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APMATH64 FUNCTION GENERATION ROUTINES 


E.1 INTRODUCTION 


This appendix presents information to help the programmer understand 
and use the function generation routines of the Advanced Math Library. 
The function generation routines are APAL64 routines that provide a 
flexible and efficient way of evaluating functions of one, two, three, 
or four variables. They do this using table lookup with linear 
interpolation. Lookup is performed by searching for the breakpoints, 
using either a binary search (successive interval halving) or a step 
search (nearest neighbor), depending on whether the user expects the 
values of the input variables to be rapidly or slowly changing from 
call to call. 


Function generation. is described in the following manner: 


Given the function F of one input variable x, for which the 
value of F is known at specific values of x (breakpoints) 
(x(1), x(2), ...), calculate the value of the function for an 
arbitrary value of x by linearly interpolating between the 
values of F at the pair of breakpoints x({i) <= x <= x(i+l). 


After determining the pair of breakpoints (x(i), x(itl)), between which 
the value of x lies, calculate the function by the following formula: 


This process is extended to two-variable functions by three 
applications of the above formula, to three-variable functions by seven 
applications, and four-variable functions by 15 applications. 


The function generation routines are listed below (refer to Appendix A 
for detailed descriptions): 


breakpoint search routines: BIN 
, STEP 


function evaluation routines: FUNI1 
FUN2 
FUN3 
FUN4 
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2 variables: + rn 4 


tw 


functions: PL(X,Y), F2(X,¥), F3(X,Y) 
3 X breakpoints: Xl, X2, X3 
4 Y breakpoints: Yl, Y2, Y3, Y4 


Coordinate value breakpoint tables: 


XBRK(1,1) = Xl YBRK(1,1) = Yl 

(2,1) = X2 (2,1) = Y2. 

(3,1) = X3 (3,1) = ¥3 

(1,2) = 1.0/(X2-X1) (4,1) = ¥4 
(2,2) = 1.0/(X3-X2) (1,2) = 1.0/(Y2-Y1) 
(3,2) = 0.0 (2,2) = 1.0/(Y3-¥2) 
(3,2) = 1.0/(Y4-Y¥3) 

Xl < X2 < X3 (4,2) = 6.6 


Yl < Y2 < ¥3 < x4 


Taken together, these two breakpoint tables specify a 3 X 4 rectangular 
grid of points in the X-Y plane. 


Function value breakpoint table: 


FBRK(1,1,1) = F1(X1,YL) 
(2,1,1) = P1L(X2,Y1) 
(3,1,1) = PL(X3,Y¥L) 
(1,2,1) = F1(X1,Y2) 
(2-2-1) = PL(X2,¥2) 
(3,2,1) = PL(X3,¥2) 
(L,3,1) = PL(X1,¥3) 
(2,3,1) = FL(X2,Y3) 
(3,3,1) = PL(X3,¥3) 
(1,4,1) = FL(X1,¥4) 
(2,4,1) = PL(X2,Y4) 
(3,4,1) = FL(X3,¥4) 
(1L,1,2) = F2(X1,%L) 
(3,4,2) = F2(X3,¥4) 
(1,1,3) = F3(X1,YL) 

: -6747- 


(3,4,3) = F3(X3,Y4) 


Figure E-1 Example Coordinate and Function Value Breakpoint Tables 
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where 


~s 
ie 
mY 
tr 
| a 
\ 


= X coordinate value of the first input point 


s 

K 

Nh 

= 
' 


= Y coordinate value of the first input point 


E.3 CALLING APMATH64 FUNCTION GENERATION ROUTINES 


The function generation package is used with System Job Executive (SJE) 
as follows: 


APFTN64 <----> Advanced Math 
driver Library routines 


The user must supply the APFTN64 driver, which contains calls to the 
appropriate Advanced Math Library routines. The coordinate value 
tables, function value table, and the input points are generated in the 
APFTNG64 driver. The APFTN64 driver routine does the following: 
@ Generates the coordinate value breakpoint tables. 
@e Generates the function value breakpoint table. 
@ Specifies the input points. 
e Sets up a loop to process the input points. 
@ For each input point, determines the appropriate breakpoint 
pair for each of the coordinates of the input point by calling 
the BIN or STEP routine for each coordinate. (This feature 


makes input point data structure arbitrary.) 


e Calls the appropriate function evaluation routine (i.e., FUNI, 
FUN2, FUN3, or FUN4 from the Advanced Math Library). 


Refer to the Advanced Math Library documentation and the individual 
program headers for descriptions of these programs. 
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The structure of the output function value array FVAL is arbitrary to 
the extent that each call to the Advanced Math Library function 
generation routine returns the interpolated values for all of the 
functions at the given input point in one array. For this reason, FVAL 
is perhaps most conveniently dimensioned FVAL(NF,NIP). 


Lines 35 through 61 load the coordinate value breakpoint tables. In 
the FUN4 example below, the program assumes the function values to be 


Lines 65 through 73 load the function value breakpoint table. In this 
example, it is done by simply cycling through all possible coordinate 
value combinations, evaluating the four functions at each point. 


Lines 77 through 198 specify the input points calling for interpolated 
values for each of the four functions. 


Lines 192 through 129 call the APMATH64 BIN and FUN4 subroutines, pass 
the tables and other arrays as arguments, and write out the results. 
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(8855) 

(8856) * WBRK(1,1)=-25.9 

(8257) WBRK(2,1)=-15.9 

(0258) WBRK(3,1)=9.9 

(8859) WBRK(1,2)=1.9/(WBRK(2,1)-WBRK(1,1)) 
(2262) WBRK(2,2)=1.9/(WBRK(3,1)-WBRK(2,1)) 
(9861) WBRK(3,2)=9.9 

(8862) 


(8863) C LOAD FBRK ARRAY 
(2864) 


(2265) DO 198 14=1,NW 


(2266) DO 199 13=1,NZ 

(0967) DO 19 I2=1,NY 

(868) DO 19 I1=1,NX 

(269) FBRK(I1,12,13,14,1)=XBRK(1I1,1)+YBRK(1I2,1)+ZBRK(I3,1)*WBRK(14,1) 
(2872) FBRK(I1,12,13,14,2)=XBRK(I1,1)*WBRK(14,1)+YBRK(I2,1)+ZBRK(13,1) 
(2871) FBRK(I1,12,13,14,3)=XBRK(I1,1)+WBRK(1I4,1)*YBRK(I2,1)+ZBRK(I3,1) 
(0272) FBRK(I1,12,13,14,4)=XBRK(I1,1)*ZBRK(I3,1)+WBRK(1I4,1)*YBRK(I2,1) 
(0973) 19@ CONTINUE 

(2274) 

(0875) C LOAD X,Y¥,Z,W ARRAYS 

(9276) 

(2277) X(1)=9.3 

(2278) ¥(1)=-5.9 

(2279) Z(1)=5.1 

(882) W(1)=-1.5 

(0481) 

(8282) X(2)=1.1 

(0883) ¥(2)=-3.8 

(2884) Z(2)=4.8 

(8885) W(2}=-2.2 

(9286) 

(2887) X(3)=0.9 

(2288) ¥(3)=-9.8 

(9889) Z(3)=7.5 

(2692) W(3)=-13.8 

(8891) 

(0292) X(4)=2.9 

(9093) ¥(4)=-6.9 

(994) Z(4)=6.0 

(2995) W(4)=-15.8 

(9996) : 

(2297) X(5)=0.4 

(8998) Y(5)=-5.9 

(2299) Z(5)=4.5 

(2192) W(5)=-7.5 

(2181) 

(2192) DO 15g I1=1,NIP 

(9193) CALL BIN(XBRK,X(I1),IX,DRX,NX) 

(194) CALL BIN(YBRK,Y(I1),IY¥,DRY,NY) 

(2195) CALL BIN(ZBRK,Z(I1),IZ,DRZ,NZ) 

(2196) CALL BIN(WBRK,W(I1),IW,DRW, NW) 

(2187) CALL FUN4(FBRK,NX,NY,NZ,NW,NF,IX,IY,1Z,IW, 
(21988) . DRX, DRY, DRZ,DRW, FVAL(1,1I1)) 
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SIMULATION LIBRARY ROUTINES 


F.1 INTRODUCTION 


The Simulation Library contains a set of routines which are useful in 
modeling various continuous systems. These continuous systems are 
characterized by ordinary differential equations (ODE) and 
three-dimensional coordinate transformations of rigid bodies, which 
Simulate physical models. 


The methods provided for solving ODE's include Runge-Kutta and Euler 
explicit methods, which require no previous evaluation of functions or 
derivatives, as well as multistep Adams implicit and explicit methods, 
which require previous evaiuation of the function and one or more 
previous derivatives. These multistep methods can be started with 
lower order methods or with the Runge-Kutta routine. Once started, the 
multistep routines require only a single evaluation of the derivative 
functions per call. The fourth order Runge-Kutta method requires four 
evaluations per time step. 


The three-dimensional rotation matrix routine forms a rotation matrix 
from a sequence of rotational specifications and can be used in 
conjunction with routine CTRN3 to perform three-dimensional coordinate 
transformations consisting of rotation plus translation. 


z 
cosine and sine of an angle, both of which are 
geometric transformations and graphic output. 


F.2 SINGLE STEP METHODS 


RKGTF Runge-Kutta-Gill-Thompson: a fourth order single step method 
to solve a system of ordinary differential equations (ODE's) 
using Thompson's numerical enhancement of the Runge-Kutta-Gill 
method. The routine requires an APFTN64 user subroutine to 
evaluate the derivatives. 


ABP1 Adams~-Bashforth predictor order one: a single step predictor 
method, also known as Euler's method, for solving ODE's. 


AMC1 Adams-Moulton corrector order one: a single step predictor 


method, also known as the backward Euler method, used for 
corrections to "stiff" ODE's. 
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APFTN64 ROUTINE FOR USE WITH RKGTF 


SUBROUTINE DFUN(T,N,Y,F) 


“Cc 
Cc *** DFUN *** SAMPLE APFTN64 ROUTINE *** 
Cc : 

DIMENSION Y(N), F(N) 
C 

DO ig I=1,N 

F(I)=Y¥(T) 
1g CONTINUE 

Cc 
Cc CORRESPONDS TO SOLUTIONS OF FORM: 
Cc 
Cc Y(I) = YO * EXP(T) 
Cc 

RETURN 

END 
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FOG RELEASE 
OLD ROUTINES 


FMMM3 2 

MMUL32 

ZVABS, VABS 
ZVADD, VADD 
ZVFLT, VFLOAT 
ZVIFIX, VIFIX 
ZVMSA, VMSA 
ZVMUL, VMUL 
ZVNEG, VNEG 
ZVRVRS, VRVRS 
ZVSADD, VSADD 
ZVSMA, VSMA 
ZVSMSA, VSMSA 
ZVSMSB, VSMSB 
ZVSMUL, VSMUL 
ZVSQ, VSQ 
ZVSUB, VSUB 
ZVSWAP, VSWAP 
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LIST OF SUPERSEDED ROUTINES 


VSADD 
VSMA 
VSMSA 
VSMSB 
VSMUL 
VsQ 
VSUB 
VSWAP 


The replacement routines for FMMM32 and MMUL32Z inciude the same 
functionality as FMMM32 and MMUL32 and are also more general. 


FS3 RELEASE 
OLD ROUTINES 


AIMAG 
CSFR 
CSFS 
CSSV 
CUFR 
CUFS 
CUSV 
EXTRU 
FLOAT 
IFIX 
INSERT 
LOC 
RSFR 
RSFS 
RSSV 
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NEW ROUTINES 


AIMAG (APFTN64 intrinsic) 


CSFR2 
CSFS2 
CSSV2 
CUFR2 
CUFS2 
CUSV2 


EXTRACT (APFTN64 intrinsic) 
FLOAT (APFTN64 intrinsic) 
IFIX (APFTN64 intrinsic) 


LOC (APFTN64 intrinsic) 


RSFR2 
RSFS2 
RSSV2 
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EXCEPTIONS ENABLED ROUTINES INFORMATION AND INTERNAL SUBROUTINES 


H.1 EXCEPTIONS ENABLED ROUTINES INFORMATION 


Beginning with the G@@ Release, all APMATH64 routines report valid 
exceptions. 


H.2 INTERNAL SUBROUTINES 


The following routines are used only as internal subroutines by other 
APMATH64 routines. These routines are listed here to facilitate 
interpretation of program tracebacks. 


INTERNAL SUBROUTINE CALLING ROUTINE(S) 

ADV2 CFFT, CFFTB, CFFTI, XCFFT 

ADV4 CFFT, CFFTB, CFFTI, XCFFT 

ALTINP CCEPS 

BITREV CFFT, CFFTI 

CBEAJY CBEJYH 

CBEDH RKGTF 

CBEDJ RKGTF 

CBERHY CBEJYH 

CBERJS CBEJYH 

CBERYH CBEJYH 

CLSTAT CFFT, CFFTB, HAMM, REALTR, STSTAT, BLKMAN, 
HANN, CFFTI, IIRELT, IREALT, XCFFT 

CTOR RFFT2D 

ENTVAR . SIMPLE 

FFT2 CFFT, CFFTB, FFT2B, CFFTI, XCFFT 

FFT2B CFFTB 

FFT4 CFFT, CFFTB, FFT4B, CFFTI 

FPT4B CFFTB 

IFFT 4 CFFTI 

IIRELT RFTII 

INTEG CBEJYH . 

IREALT RFFTI, RFTII, IIRELT 

LPSPFI SIMPLE 

PHAUNW CCEPS 

PHCHCK PHAUNW 

REALTR RFFT, RFFTB, RFFTI, RFTII 

RKGTF INTEG 

RTOC RFFT2D 

SET2 4B FFT2B, FFT4B 

SPCVAL PHAUNW 

STSTAT CFFT, CFFTB, HAMM, REALTR, STSTAT, BLKMAN, 
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CCMMUL 
CDET 
CDOTPR 
CFFT 
CFFTB 
CrFTM 
CFFTSC 
CGMMUL 
CMATIN 
CMDET 
CMFACT 
CMMTRC 
CMMUL 
CMSOLV 
CMTRAC 
CMTRAN 
CMVML3 
CMVML4 
CONV 


CDMMTITT. 


ee ee ed 


CROSSP 
CRVADD 
CRYDIV 
CRVMUL 
CRVSUB 
CSOLV 
CSOLVQ 
CTRN2 
CTRN3 
CVABS 
CVADD 
CVCOMB 
CVCONJ 
CVEXP 
CVFILL 
CVMA 
CVMAGS 
CVMEXP 
CYVMOV 
CVMUL 
CVNEG 
CVRCIP 
CVREAL 


APMATH64 ROUTINES IN PAGE ORDER AND BY TYPE 


COMPLEX 
COMPLEX 
COMPLEX 
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BASIC MATH LIBRARY (VOLUME 1) 


MATRIX MULTIPLY 
MATRIX DETERMINANT 
DOT PRODUCT 


COMPLEX-TO-COMPLEX FFT (IN PLACE) 
COMPLEX~TO-COMPLEX FFT (NOT IN PLACE) 
MIXED-RADIX COMPLEX FFT (NOT-IN-PLACE) 


COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 


FFT SCALE 

GENERAL MATRIX MULTIP 
MATRIX INVERSE 

MATRIX DETERMINANT 
MATRIX L/U FACTORIZATION 
MATRIX MULTIPLY TRACE 
MATRIX MULTIPLY 

MATRIX EQUATION SOLVER 
SUB-MATRIX TRACE 

SUB-MATRIX TRANSPOSE 

3X3 MATRIX MULT. 3D VECTORS 
4X4 MATRIX MULT. 4D VECTORS 


t 


CONVOLUTION (CORRELATION) 


COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 


3D CROSS PRODUCT 

AND REAL VECTOR ADD 

AND REAL VECTOR DIVIDE 
AND REAL VECTOR MULTIPLY 
AND REAL VECTOR SUBTRACT 
SYSTEM SOLVER 

MATRIX EQUATION SOLVER 


2-D COORDINATE TRANSFORM 
3-DIMENSIONAL COORDINATE TRANSFORMATION 


COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 


VECTOR ABSOLUTE VALUE 
VECTOR ADD 

VECTOR COMBINE 

VECTOR CONJUGATE 

VECTOR EXPONENTIAL 
VECTOR FILL 

VECTOR MULTIPLY AND ADD 
VECTOR MAGNITUDE SQUARED 
VECTOR MULTIPLY EXPONENTIAL 
VECTOR MOVE 

VECTOR MULTIPLY 

VECTOR NEGATE 

VECTOR RECIPROCAL 


FORM COMPLEX VECTOR OF REALS 
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BASIC MATH LIBRARY (cont. ) 
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RMSQV 
SCJMA 
SGEFA 
SGESL 
SGTSL 
SMMM 
SMMMV 
SN2 
SOLVEQ 
STMM 
SVE 
SVEMG 
SVESQ 
SVS 
TRIDIA 
VAAM 
VABS 
VACOS 
VADD 
VAINT 
VALG 
VALG1Z 
VAM 
VASIN 
VASM 
VATAN: 


tT7AMARTD 
VEIL FAIS 


VCLIP 
VCLR 
vcos 
VCOSH 
VDIV 
VEUCL2 
VEXP 
VEXP19 
VFILL 
VFRAC 
VIABS 
VIADD 
VICLIP 
VIDIV 
VIMAG 
VIMUL 
VINDEX 
VINEG 


od et er a 


ROOT-MEAN-SQUARE OF VECTOR ELEMENTS 
SELF-CONJUGATE MULTIPLY AND ADD 
REAL -GENERAL MATRIX FACTOR 

REAL GENERAL MATRIX SOLVE 


TRIDIAGONAL MATRIX SOLVER 
SUBMATRIX MULTIPLY 


dee EE 2 oe io hcl oD ae ee ee ote 


SUBMATRIX MULTIPLY 
SQUARED DISTANCE BETWEEN TWO VECTORS 


LINEAR 


EQUATION SOLVER 


SUBMATRIX TRANSPOSE & MULTIPLY 


SUM OF 
SUM OF 
SUM OF 
SUM OF 


VECTOR ELEMENTS 

VECTOR ELEMENT MAGNITUDES 
VECTOR ELEMENT SQUARES 
VECTOR SIGNED SQUARES 


~~. mses = 2 RAS ewer 


TRIDIAGONAL MATRIX SOLVER 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


VEMMNAD 
V bee BEN 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


ADD, ADD, AND MULTIPLY 
ABSOLUTE VALUE 
ARCCOSINE 

ADD 

TRUNCATE 

LOGARI THM 

BASE 19 LOGARITHM 

ADD AND MULTIPLY 
ARCSINE 

ADD AND SCALAR MULTIPLY 


ARCTANGENT 
ARCTANGENT (2 ARGUMENTS) 


cadre TNS deta ose We steht wD 


CLIP 

CLEAR 

COSINE 

COSINE (HYPERBOLIC) 
DIVIDE 

EUCLIDEAN DISTANCE 
EXPONENTIAL 
EXPONENTIAL (BASE 19) 
FILL 

TRUNCATE TO FRACTION 


ABSOLUTE VALUE 
INTEGER ADD 
INVERTED CLIP 
INTEGER DIVIDE 


EXTRACT IMAGINARIES OF COMPLEX VECTOR 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


INTEGER MULTIPLY 
INDEX 

INTEGER NEGATE 
INTEGER SUBTRACT 
LOGICAL ADD 

LOGICAL EQUIVALENCE 
LIMIT 


LOGICAL VECTOR MERGE 
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APPENDIX I 


BASIC MATH LIBRARY (cont.) 


VTSADD VECTOR TM SCALAR ADD A - 239 
VTSMA VECTOR TM SCALAR MULTIPLY AND ADD A - 231 
VTSMUL VECTOR TM SCALAR MULTIPLY A - 232 
ADVANCED MATH LIBRARY (VOLUME 2) 
CH COMPLEX HERMITIAN EIGENSYSTEM SOLVER A 234 
EIGRS REAL SYMMETRIC EIGENSYSTEM SOLVER A 237 
HTRIBK COMPLEX HERMITIAN EIGENVECTORS A 239 
HTRIDI COMPLEX HERMITIAN TRIDIAGONALIZATION A 241 
IMTQL1 DIAGONALIZE TRIDIAGONAL MATRIX A 243 
IMTQL2 DIAGONALIZE A TRIDIAGONAL MATRIX A 245 
RS REAL SYMMETRIC EIGENSYSTEM SOLVER A 247 
SIMPLE REVISED SIMPLEX A 249 
SKYSOL SKYLINE FORMAT EQUATION SOLVER A 254 
TRED1 TRIDIAGONALIZE SYMMETRIC MATRX A 256 
TRED2 TRIDIAGONALIZE A SYMMETRIC MATRIX A 258 
VASORT VECTOR SORT ALGEBRAIC VALUES A 259 
VISORT VECTOR SORT INTEGER VALUES A 269 
VSORT VECTOR SORT WITH INDICES A 261 
SIGNAL PROCESSING LIBRARY 
ACORF AUTO-CORRELATION (FREQUENCY-DOMAIN) A 264 
ACORT AUTO-CORRELATION (TIME-DOMAIN) A 265 
ASPEC ACCUMULATING AUTO-SPECTRUM A 266 
BLKMAN BLACKMAN WINDOW MULTIPLY A 267 
CCEPS PHASE UNWRAP AND COMPLEX CEPSTRUM A 268 
CCORF CROSS~-CORRELATION (FREQUENCY-DOMAIN) A 272 
CCORT CROSS~-CORRELATION (TIME-DOMAIN) A 273 
CFFTI COMPLEX FFT WITH INTERPOLATION A 274 
COHER COHERENCE FUNCTION A 275 
CSPEC ACCUMULATING CROSS-SPECTRUM A 276 
DECFIR DECIMATION A 277 
ENVEL ENVELOPE DETECTOR A 279 
HAMM HAMMING WINDOW MULTIPLY A 288 
HANN HANNING WINDOW MULTIPLY A 281 
HIST HISTOGRAM A 282 
HLBRT HILBERT TRANSFORMER A 283 
LPAUTO LINEAR PREDICTION AUTOCORRELATION A 284 
PKVAL PEAK AND VALLEY PICKING A 286 
RDFT REAL DISCRETE FOURIER TRANSFORM A 288 
RFFTI REAL FFT WITH INTERPOLATION A 299 
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LINPACK BLAS LIBRARY (cont.) 


SROT PLANE ROTATION A - 355 
SROTG GIVENS PLANE ROTATION A - 356 
SROTM MODIFIED GIVENS ROTATIONS A - 357 
SROTMG MODIFIED- GIVENS PLANE ROTATION SETUP A - 358 
SSCAL REAL SCALAR TIMES VECTOR A - 369 
SSWAP INTERCHANGES VECTORS A - 361 
SIMULATION LIBRARY 
ABP1 ADAMS-BASHFORTH PREDICTOR (ORDER 1) A - 363 
ABP2 ADAMS-BASHFORTH PREDICTOR (ORDER 2) A - 364 
ABP3 ADAMS-BASHFORTH PREDICTOR (ORDER 3) A - 365 
ABP4 ADAMS-BASHFORTH PREDICTOR (ORDER 4) A - 366 
ADAMS 4 ADAMS VARIABLE STEP INTEG.(ORD 4) A - 368 
AMC1 ADAMS~-MOULTON CORRECTOR (ORDER 1) A - 371 
AMC2 ADAMS~-MOULTON CORRECTOR (ORDER 2) A - 372 
AMC3 ADAMS-MOULTON CORRECTOR (ORDER 3) A - 373 
AMC 4 ADAMS-MOULTON CORRECTOR (ORDER 4) A - 374 
BIN BINARY SEARCH A - 376 
CBEIK COMPLEX BESSEL I AND K A - 378 
CBEJYH COMPLEX BESSEL J, Y, AND H A - 381 
FUN FUNCTION OF ONE VARIABLE A - 384 
FUN2 FUNCTION OF TWO VARIABLES A - 386 
FUN3 _ FUNCTION OF THREE VARIABLES A - 388 
FUN4 FUNCTION OF FOUR VARIABLES A - 399 
RKGIL RUNGE~-KUTTA-GILL INTEGRATION A - 392 
RKGTF R-K-GILL-THOMPSON INTEG. (ORDER 4) A - 395 
ROT3 3D ROTATION MATRIX, 3-ANGLE A - 397 
scsl SCALAR COS/SIN, TM INTERP.(ORD 1) A - 399 
STEP STEP SEARCH A - 499 


GEOPHYSICAL LIBRARY 


CONNMO NMO WITH CONSTANT VELOCITY A - 493 
IIR3G RECURSIVE FILTER A - 485 
KSMLV K-TH SMALLEST ELEMENT IN VECTOR A - 4987 
NMOLI NMO LINEAR INTERPOLATION A - 468 
NMOQI NMO QUADRATIC INTERPOLATION A - 419 
POST64 POST BITS TO RASTER A - 412 
RESNMO RESIDUAL NORMAL MOVEOUT A - 413 
TMCONV CONVOLUTION (CORRELATION } A - 415 
Vg1 VECTOR ZERO TRENDS A - 417 
VARNMO NMO WITH VARIABLE VELOCITY A - 418 
VRNAVG VECTOR RUNNING AVERAGE A - 429 
VSCAND VECTOR SCAN FOR ZEROS A - 421 


~] 
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TABLE MEMORY LIBRARY (cont. ) 


TTTSUB VECTOR SUBTRACT (TM-TM TO TM) A - 599 
TTVLC2 VECTOR LINEAR COMBINATION A - 591 
TVCLR TABLE MEMORY VECTOR CLEAR A - 592 
SPECIAL UTILITIES LIBRAR 
EXTRS EXTRACT A SIGNED BIT-FIELD A - 594 
PEEK MEMORY FETCH A ~ 595 
POKE STORE INTO MEMORY A - 596 
DATA FORMATTING LIBRARY 
VFLOAT CONVERT INTEGER TO FLOATING-POINT A - 598 
VIFIX VECTOR INTEGER FIX A - 599 
VPK8 VECTOR 8-BIT BYTE PACK A - 519 
VPK16 VECTOR 16-BIT BYTE PACK A - 511 
VPK32 VECTOR 32-BIT BYTE PACK A - 512 
VPKI32 VECTOR 32-BIT INTEGER PACK A - 513 
VPKR32 VECTOR REAL HALFWORD PACK A - 514 
VSCALE VECTOR SCALE AND FIX A - 515 
VSCSCL VECTOR SCAN SCALE AND FIX A - 516 
VSHFX VECTOR SHIFT AND FIX A - 517 
VSMAFX VECTOR SCALAR MULTIPLY, ADD AND FIX A - 518 
VUP8 VECTOR 8-BIT BYTE UNPACK A - 519 
VUP16 VECTOR 16-BIT BYTE UNPACK A - 529 
VUP32 VECTOR 32-BIT BYTE UNPACK A - 521 
VUPR32 VECTOR HALFWORD REAL UNPACK A - 522 
VUPS8 VECTOR 8-BIT SIGNED BYTE UNPACK A - 523 
VUPS16 VECTOR 16-BIT SIGNED BYTE UNPACK A - 524 
VUPS32 VECTOR 32-BIT SIGNED BYTE UNPACK A - 525 
VUSI32 VECTOR 32-BIT SIGNED INTEGER UNPACK A - 526 
VUUI32 VECTOR 32-BIT UNSIGNED UNPACK A - 527 
DOUBLE PRECISION LIBRARY 
DADD DOUBLE TO DOUBLE-PRECISION ADD A - 529 
DADDRR SINGLE TO DOUBLE-PRECISION ADD A - 538 
DADOT DOUBLE ACCUMULATE DOT PRODUCT A - 531 
DDOTRR DOUBLE DOT PRODUCT REAL REAL A - 532 
DMUL DOUBLE TO DOUBLE-PRECISION MULTIPLY A - 533 
DMULRR SINGLE TO DOUBLE PRECISION MULTIPLY A - 534 
DNEG NEGATE DOUBLE-PRECISION NUMBER A - 535 
DSUB DOUBLE TO DOUBLE-PRECISION SUBTRACT A - 536 
DSUBRR SINGLE TO DOUBLE-PRECISION SUBTRACT A - 537 
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APPENDIX J 


APMATH64 ROUTINES IN ALPHABETICAL ORDER 


NAME DESCRIPTION PAGE 

ABP1 ADAMS~BASHFORTH PREDICTOR (ORDER 1) A - 363 
ABP2 ADAMS~-BASHFORTH PREDICTOR (ORDER 2) A - 364 
ABP3 ADAMS~-BASHFORTH PREDICTOR (ORDER 3) A - 365 
ABP4 ADAMS~BASHFORTH PREDICTOR (ORDER 4) A - 366 
ABS REAL NUMBER ABSOLUTE VALUE A - 539 
ACORF AUTO-CORRELATION (FREQUENCY-DOMAIN) A - 264 
ACORT AUTO-CORRELATION (TIME-DOMAIN) A - 265 
ACOS REAL NUMBER ARCCOSINE A - S42 
ADAMS4 ADAMS VARIABLE STEP INTEG.(ORD 4) A - 368 
AINT TRUNCATE REAL NUMBER A - 541 
ALOG19 REAL NUMBER LOGARITHM A - 543 
ALOG REAL NUMBER LOGARITHM A - 542 
AMC1 ADAMS~MOULTON CORRECTOR (ORDER 1) A - 371 
AMC2 ADAMS~MOULTON CORRECTOR (ORDER 2) A - 372 
AMC3 ADAMS~-MOULTON CORRECTOR (ORDER 3) A - 373 
AMC4 ADAMS-MOULTON CORRECTOR (ORDER 4) A - 374 
AMOD REAL/REAL DIVIDE REMAINDER A - 544 
ANINT ROUND REAL NUMBER TO NEAREST WHOLE A - 545 
ASIN REAL NUMBER ARCSINE A - 546 
ASPEC ACCUMULATING AUTO-SPECTRUM A - 266 
ATAN2 ARCTANGENT OF RATIO OF REAL NUMBERS A - 548 
ATAN ARCTANGENT OF REAL NUMBER A - 547 
BIN BINARY SEARCH A - 376 
BLKMAN BLACKMAN WINDOW MULTIPLY A - 267 
CABS COMPLEX NUMBER ABSOLUTE VALUE A - 549 
CAXPY COMPLEX A * X + ¥ A - 326 
CAXPYN NESTED COMPLEX A * X + Y A - 327 
CBEIK COMPLEX BESSEL I AND K A - 378 
CBEJYH COMPLEX BESSEL J, Y, AND H A - 381 
CCEPS PHASE UNWRAP AND COMPLEX CEPSTRUM A - 268 
CCMMUL COMPLEX MATRIX MULTIPLY A- 2 

CCOPY COMPLEX VECTOR COPY A - 339 
CCORF CROSS-CORRELATION (FREQUENCY-DOMAIN) A - 272 
CCORT CROSS~-CORRELATION (TIME-DOMAIN) A - 273 
ccos COMPLEX NUMBER COSINE A - 559 
CDET COMPLEX MATRIX DETERMINANT A- 4 

CDIV COMPLEX/COMPLEX DIVIDE A - 551 
CDIVCR REAL/COMPLEX DIVIDE A - 552 
CDIVRC COMPLEX/REAL DIVIDE A - 553 
CDOTC COMPLEX INNER PRODUCT A - 331 
CDOTN NESTED COMPLEX DOT PRODUCT A - 332 
CDOTPR COMPLEX DOT PRODUCT A- 6 

CDOTU COMPLEX DOT PRODUCT A - 335 
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CUFS2 
cusv2 
CVABS 
CVADD 
CVCOMB 
CVCONT 
CVEXP 
CVFILL 
CVMA 
CVMAGS 
CVMEXP 
CVMOV 
CVMUL 
CVNEG 
CVRCIP 
CVREAL 
CVSMA 
CVSMUL 
CVSUB 
DADD 
DADDRR 
DADOT 
DDOTRR 
DECFIR 
DEQ22 
DIM 
DMUL 
DMULRR 
DNEG 
DOTPR 
DSUB 
DSUBRR 
EIGRS 


FUNL 
FUN2 
FUN3 
FUN4 
GENTAB 
GRAD2D 
GRD2DB 


SPARSE COMPLEX UNSYM FACTOR & SOLVE 
SPARSE COMPLEX UNSYMMETRIC SOLVE 
COMPLEX VECTOR ABSOLUTE VALUE 
COMPLEX VECTOR ADD 

COMPLEX VECTOR COMBINE 

COMPLEX VECTOR CONJUGATE 

COMPLEX VECTOR EXPONENTIAL 

COMPLEX VECTOR FILL 

COMPLEX VECTOR MULTIPLY AND ADD 
COMPLEX VECTOR MAGNITUDE SQUARED 
COMPLEX VECTOR MULTIPLY EXPONENTIAL 
COMPLEX VECTOR MOVE 

COMPLEX VECTOR MULTIPLY 

COMPLEX VECTOR NEGATE 

COMPLEX VECTOR RECIPROCAL 

FORM COMPLEX VECTOR OF REALS 
COMPLEX VECTOR SCALAR MULTIPLY AND ADD 
COMPLEX VECTOR SCALAR MULTIPLY 
COMPLEX VECTOR SUBTRACT 

DOUBLE TO DOUBLE-PRECISION ADD 
SINGLE TO DOUBLE-PRECISION ADD 
DOUBLE ACCUMULATE DOT PRODUCT 
DOUBLE DOT PRODUCT REAL REAL 
DECIMATION 

DIFFERENCE EQUATION, 2 POLES, 2 ZEROS 
REAL/REAL POSITIVE DIFFERENCE 
DOUBLE TO DOUBLE-PRECISION MULTIPLY 
SINGLE TO DOUBLE PRECISION MULTIPLY 
NEGATE DOUBLE-PRECISION NUMBER 

DOT PRODUCT 

DOUBLE TO DOUBLE-PRECISION SUBTRACT 
SINGLE TO DOUBLE-PRECISION SUBTRACT 
REAL SYMMETRIC EIGENSYSTEM SOLVER 
ENVELOPE DETECTOR 

EXPONENTIAL OF REAL NUMBER 

EXTRACT A SIGNED BIT-FIELD 

FAST MATRIX MULTIPLY 

FAST MATRIX MULTIPLY 

FUNCTION OF ONE VARIABLE 

FUNCTION OF TWO VARIABLES 

FUNCTION OF THREE VARIABLES 
FUNCTION OF FOUR VARIABLES 

GENERATE TWIDDLE FACTOR TABLE 
MAXIMUM GRADIENT FILTER 

MAXIMUM GRADIENT FILTER WITH BOUND 
HAMMING WINDOW MULTIPLY 

HANNING WINDOW MULTIPLY 

HISTOGRAM 

HILBERT TRANSFORMER 

COMPLEX HERMITIAN EIGENVECTORS 
COMPLEX HERMITIAN TRIDIAGONALIZATION 
INTEGER ABSOLUTE VALUE 

INDEX OF LARGEST COMPLEX ELEMENT 
INTEGER/INTEGER POSITIVE DIFFERENCE 
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PAS3F 
PAS3I 
PAS4P 
PAS4I 
PASSF 
PASSI 
PEEK 
PFINV 
PKVAL 
POKE 
POLAR 
POST64 
RAN 
RCMMUL 


RADIX-3 
RADIX-3 
RADIX~4 
RADIX-4 
RADIX-5 
RADIX~-5 


FORWARD 
INVERSE 
FORWARD 
INVERSE 
FORWARD 
INVERSE 


COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 
COMPLEX 


FFT 
FPT 
FFT 
FFT 
FFT 
FFT 


PASS 
PASS 
PASS 
PASS 
PASS 
PASS 


MEMORY FETCH 

MATRIX INVERSE (PRODUCT FORM) 

PEAK AND VALLEY PICKING 

STORE INTO MEMORY 

RECTANGULAR TO POLAR CONVERSION 
POST BITS TO RASTER . 

SCALAR RANDOM NUMBER GENERATOR 
REAL-COMPLEX MATRIX MULTIPLY 

REAL DISCRETE FOURIER TRANSFORM 
REAL/REAL DIVIDE 

POLAR TO RECTANGULAR CONVERSION 
RESIDUAL NORMAL MOVEOUT 

REAL TO COMPLEX 2-DIMENSIONAL FFT 
REAL-TO-COMPLEX FFT (IN PLACE) 
REAL-TO-COMPLEX FFT (NOT IN PLACE) 
REAL FFT WITH INTERPOLATION 
MIXED-RADIX REAL FFT (NOT-IN-PLACE) 
REAL FFT SCALE AND FORMAT 

REAL FFT WITH QUARTER INTERPOLATION 
REAL GENERAL MATRIX MULTIPLY 
RUNGE-KUTTA-GILL INTEGRATION 
R-K-GILL-THOMPSON INTEG. (ORDER 4) 
ROOT-MEAN-SQUARE OF VECTOR ELEMENTS 
3D ROTATION MATRIX, 3-ANGLE 

REAL TO REAL POWER 

REAL TO INTEGER POWER 

REAL RECIPROCAL 

REAL SYMMETRIC EIGENSYSTEM SOLVER 
SPARSE REAL SYMMETRIC FACTOR 
SPARSE REAL SYMM FACTOR & SOLVE 
RECIPROCAL SQUARE ROOT 

SPARSE REAL SYMMETRIC SOLVE 
SPARSE REAL UNSYMMETRIC FACTOR 
SPARSE REAL UNSYM FACTOR & SOLVE 
SPARSE REAL UNSYMMETRIC SOLVE 

SUM OF MAGNITUDES 

REAL A * X + Y 

NESTED REAL A * X + Y 

SUM OF REAL AND IMAGINARY MAGNITUDES 
SELF-CONJUGATE MULTIPLY AND ADD 
COMPLEX EUCLIDEAN NORM 

VECTOR COPY 

SCALAR COS/SIN, TM INTERP.(ORD 1) 
DOT PRODUCT OF REAL VECTORS 
NESTED REAL DOT PRODUCT 

SPARSE DOT PRODUCT 

REAL GENERAL MATRIX FACTOR 

REAL GENERAL MATRIX SOLVE 
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TTTSUB 
TTVLC2 
TVCLR 
VG1 
VAAM 
VABS 
VACOS 
VADD 
VAINT 
VALG1S 
VALG 
VAM 
VARNMO 
VASIN 
VASM 
VASORT 
VATAN2 
VATAN 
VAVEXP 
VAVLIN 
VCLIP 
VCLR 
VCcos 
VCOSH 
VDBPWR 
VDIV 
VEUCL2 
VEXP19 
VEXP 
VFILL 
VFLOAT 
VFRAC 
VIABS 
VIADD 
VICLIP 
VIDIV 
VIFIX 
VIMAG 
VIMUL 
VINDEX 
VINEG 
VISORT 
VISUB 
VLAND 
VLEQV 
VLIM 


VECTOR 
VECTOR 


SUBTRACT (TM-TM TO TM) 
LINEAR COMBINATION 


TABLE MEMORY VECTOR CLEAR 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


ZERO TRENDS 

ADD, ADD, AND MULTIPLY 
ABSOLUTE VALUE 
ARCCOSINE 

ADD 

TRUNCATE 

BASE 19 LOGARITHHM 
LOGARITHM 

ADD AND MULTIPLY 


NMO WITH VARIABLE VELOCITY 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


CONVERT INTEGER TO FLOATING-POI 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


ARCSINE 

ADD AND SCALAR MULTIPLY 
SORT ALGEBRAIC VALUES 
ARCTANGENT (2 ARGUMENTS) 
ARCTANGENT 

EXPONENTIAL AVERAGING 
LINEAR AVERAGING 

CLIP 

CLEAR 

COSINE 

COSINE (HYPERBOLIC) 
CONVERSION TO DB (POWER) 
DIVIDE 

EUCLIDEAN DISTANCE 
EXPONENTIAL (BASE 12) 
EXPONENTIAL 

FILL 

TRUNCATE TO FRACTION 
ABSOLUTE VALUE 

INTEGER ADD 

INVERTED CLIP 

INTEGER DIVIDE 

INTEGER FIX 


EXTRACT IMAGINARIES OF COMPLEX VECTOR 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


INTEGER MULTIPLY 
INDEX 

INTEGER NEGATE 

SORT INTEGER VALUES 
INTEGER SUBTRACT 
LOGICAL ADD 

LOGICAL EQUIVALENCE 
LIMIT 


LOGICAL VECTOR MERGE 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


VECTOR 
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LOGICAL OR 

LOGICAL EXCLUSIVE GOR 
MULTIPLY AND ADD 
MAXIMUM 

MAXIMUM MAGNITUDE 
MINIMUM 
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VTSMA 
VTSMUL 
VUP16 
VUP32 
VUP8 
VUPR32 
VUPS16 
VUPS32 
VUPS8 
VUST32Z 
VUUI32 
VXCS 
WIENER 


VECTOR 
VECTOR 
VECTOR 
VECTOR 
VECTOR 


VECTOR 
VECTOR 
WIENER 
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TM SCALAR MULTIPLY AND ADD 
TM SCALAR MULTIPLY 

16-BIT BYTE UNPACK 

32-BIT BYTE UNPACK 

8-BIT BYTE UNPACK 

HALFWORD REAL UNPACK 
16-BIT SIGNED BYTE UNPACK 
32-BIT SIGNED BYTE UNPACK 
8-BIT SIGNED BYTE UNPACK 
32-BIT SIGNED INTEGER UNPACK 
32-BIT UNSIGNED UNPACK 
MULTIPLIED BY SIN AND COS 
LEVINSON ALGORITHM 
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APMATH64 KEY WORD INDEX 


index of APMATH64 routines is sorted by key words that appear in 
routine title. Each title can contain more than one key word. 

The key words are listed alphabetically to the right of the gap running 
down the center of each page. 


To use the key word index, locate a key word that is representative of 
the desired APMATH64 function. Applicable APMATH64 routine names and 
titles can be found on the same line with each occurrence of the key 
word. The routine name appears in brackets ({ ]). The routine title 
immediately follows the routine name and continues on the other side of 
the gap when necessary. The ellipsis (...}) is placed directly after 
the last word in the title if the line wraps around. The page where a 


particular routine is documented can be found in Appendix J. 


[VSIMPS]VECTOR SIMPSON'S 
EXPONENTIAL BASE 
[VALG1Z]VECTOR BASE 

[ VPK16 ] VECTOR 

{ VUP16 ] VECTOR 

[ VUPS16 ] VECTOR 

[ CSROT ]JCOMPLEX 

[CFFT2D ]COMPLEX TO COMPLEX 
(RFFT2D]REAL TO COMPLEX 
(ROT3]3D ROTATION MATRIX, 
(VPK32 ] VECTOR 

( VUP32 ] VECTOR 

[ VPKI32 ]VECTOR 

[ VUPS32 ] VECTOR 
[VUSI32 ] VECTOR 

{ VUUI32 ]VECTOR 

[CROSSP ]COMPLEX 

3X3 MATRIX MULT. 

[MVML3 ]MATRIX VECTOR MULTIPLY 
[ CMVML3 ] COMPLEX 

4X4 MATRIX MULT. 

(MVML4 ]MATRIX VECTOR MULTIPLY 
{ CMVML4 ] COMPLEX 

[ VPK8 ] VECTOR 

[ VUP8 ] VECTOR 

[ VUPS8 ] VECTOR 

[ABS]REAL NUMBER 

{CABS ]JCOMPLEX NUMBER 
[CVABS ]COMPLEX VECTOR 

{ [ABS ] INTEGER 

[ ISAMAX]INDEX OF MAXIMUM 
[ VABS ] VECTOR 

[ VIABS ] VECTOR 

[ DADOT ] DOUBLE 
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1/3 RULE INTEGRATION 

1g... [VEXP19]VECTOR 

18 LOGARITHM 

16-BIT BYTE PACK 

16-BIT BYTE UNPACK 

16-BIT SIGNED BYTE UNPACK 
2-D ROTATION 

2-DIMENSIONAL FFT 
2~DIMENSIONAL FFT 

3-ANGLE 

32-BIT BYTE PACK 

32-BIT BYTE UNPACK 

32-BIT INTEGER PACK 

32-BIT SIGNED BYTE UNPACK 
32-BIT SIGNED INTEGER UNPACK 
32-BIT UNSIGNED UNPACK 

3D CROSS PRODUCT 

3D VECTORS. ..{CMVML3 ] COMPLEX 
3X3 

3X3 MATRIX MULT. 3D VECTORS 
4D VECTORS. ..[{CMVML4 ] COMPLEX 
4X4 

4X4 MATRIX MULT. 
8-BIT BYTE PACK 
8-BIT BYTE UNPACK 
8-BIT SIGNED BYTE UNPACK 
ABSOLUTE VALUE 

ABSOLUTE VALUE 

ABSOLUTE VALUE 

ABSOLUTE VALUE 

ABSOLUTE VALUE 

ABSOLUTE VALUE 

ABSOLUTE VALUE 
ACCUMULATE DOT PRODUCT 


4D VECTORS 


INDEX 


[VUP16]VECTOR 16-BIT 
{VUP32]VECTOR 32-BIT 
{VUP8]VECTOR 8-BIT 
(VUPS16]VECTOR 16-BIT SIGNED 
[VUPS32]VECTOR 32-BIT SIGNED 
(VUPS8]VECTOR 8-BIT SIGNED 
UNWRAP AND COMPLEX 
(TVCLR]TABLE MEMORY VECTOR 
[ VCLR ] VECTOR 

[ VCLIP ] VECTOR 
(VICLIP]VECTOR INVERTED 
[TMVLC2]VECTOR LINEAR 
[TTVLC2]VECTOR LINEAR 

{ CVCOMB ]COMPLEX VECTOR 
[CFFT2D]COMPLEX TO 
(RFFT2D]REAL TO 

{[CAXPYN ]NESTED 
[CCEPS]PHASE UNWRAP AND 

( CDOTN ] NESTED 

[ ICAMAX ] INDEX OF LARGEST 

{ CFFTM ]MIXED-RADIX 
[PAS2F]RADIX-2 FORWARD 
{PAS2I]RADIX-2 INVERSE 
(PAS3F ]RADIX-3 FORWARD 

{ PAS3I ]RADIX-3 INVERSE 
[PAS4F]RADIX-4 FORWARD 
{PAS4I ]RADIX-4 INVERSE 
{PASSF]RADIX-5 FORWARD 
{PASSI]RADIX-5 INVERSE 
[CEXP JEXPONENTIAL OF 
[CONJG ]CONJUGATE OF 
[CSQRT]SQUARE ROOT OF 
[CPOWRC]REAL TO 
(CPOW]COMPLEX TO 
(CSFS2 ] SPARSE 

[CSFR2 JSPARSE 

(CSSV2 ]SPARSE 

[ CUFS2 ]SPARSE 
[CUFR2 ] SPARSE 

[CUSV2 ]SPARSE 
(VIMAG]EXTRACT IMAGINARIES OF 
{VREAL]EXTRACT REALS OF 

{ CVREAL ] FORM 

(CSSCAL]REAL TIMES 

[CVCONJ ]JCOMPLEX VECTOR 
[CONNMO ]NMO WITH 

[POLAR ]RECTANGULAR TO POLAR 
{RECT]POLAR TO RECTANGULAR 
[ VDBPWR ] VECTOR 

{CONV2D j2-D 

[ TCONV ] POST-TAPERED 

[CTRN2 ]2-D 

{ CTRN3 ] 3-DIMENSIONAL 
{CCOPY ]COMPLEX VECTOR 
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BYTE UNPACK 
BYTE UNPACK 

BYTE UNPACK 

BYTE UNPACK 

BYTE UNPACK 

BYTE UNPACK 

CEPSTRUM. ..[CCEPS ] PHASE 
CLEAR 

CLEAR 

CLIP 

CLIP 

COMBINATION 

COMBINATION 

COMBINE 

COMPLEX 2-DIMENSIONAL FFT 
COMPLEX 2-DIMENSIONAL FFT 
COMPLEX A * X + Y¥ 

COMPLEX CEPSTRUM 

COMPLEX DOT PRODUCT 
COMPLEX ELEMENT 

COMPLEX FFT NOT-IN-PLACE 
COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX FFT PASS 

COMPLEX NUMBER 

COMPLEX NUMBER 

COMPLEX NUMBER 

COMPLEX POWER 

COMPLEX POWER 

COMPLEX SYMM FACTOR & SOLVE 
COMPLEX SYMMETRIC FACTOR 
COMPLEX SYMMETRIC SOLVE 
COMPLEX UNSYM FACTOR & SOLVE 
COMPLEX UNSYMMETRIC FACTOR 
COMPLEX UNSYMMETRIC SOLVE 
COMPLEX VECTOR 

COMPLEX VECTOR 

COMPLEX VECTOR OF REALS 
COMPLEXES 

CONJUGATE 

CONSTANT VELOCITY 
CONVERSION 

CONVERSION 

CONVERSION TO DB POWER 
CONVOLUTION AND CORRELATION 
CONVOLUTION CORRELATION 
COORDINATE TRANSFORM 
COORDINATE TRANSFORMATION 
COPY 


INDEX 


[ DNEG ] NEGATE 
[DSUBRR]SINGLE TO 
[DSUB]DOUBLE TO 
(CH]COMPLEX HERMITIAN 
{EIGRS]REAL SYMMETRIC 
{RS]REAL SYMMETRIC 
[HTRIBK]COMPLEX HERMITIAN 
OF LARGEST COMPLEX 
(KSMLV]K-TH SMALLEST 
“ert 
J 
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[MAAMGV JMAXLMUM MAGNITUDE 


{ MAXV ]MAXIMUM 

[MINMGV ]MINIMUM MAGNITUDE 
{MINV ]MINIMUM 

[MEAMGV ]JMEAN OF VECTOR 
[SVEMG]SUM OF VECTOR 
{MEASQV ]MEAN OF VECTOR 
[SVESQ]SUM OF VECTOR 
[MEANV ]MEAN VALUE OF VECTOR 
OF VECTOR 
{SVE]SUM OF VECTOR 
[VSUM]VECTOR SUM OF 
(LVEQ]LOGICAL VECTOR 
VECTOR GREATER THAN OR 
[LVNE]LOGICAL VECTOR NOT 
{ CMSOLV ]COMPLEX MATRIX 
{CSOLVQ]COMPLEX MATRIX 
(SKYSOL]SKYLINE FORMAT 

{ SOLVEQ ] LINEAR 

[ DEQ22 ] DIFFERENCE 
(VLEQV]VECTOR LOGICAL 
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LVOoUULIS J VOULUN 


{ SCNRM2 ] COMPLEX 
[VPOLY]VECTOR POLYNOMIAL 
{[VLXOR]VECTOR LOGICAL 

(CVEXP ]COMPLEX VECTOR 
VECTOR MULTIPLY 

{ VEXP ] VECTOR 

{ VAVEXP ] VECTOR 

{ VEXP19 ] VECTOR 

COMPLEX SYMMETRIC 

COMPLEX UNSYMMETRIC 
[RSFR2]SPARSE REAL SYMMETRIC 
[RUFR2]SPARSE REAL UNSYMMETRIC 
[SGEPA]REAL GENERAL MATRIX 
(CSFS2]SPARSE COMPLEX SYMM 
{CUFS2]SPARSE COMPLEX UNSYM 
[RSFS2]SPARSE REAL SYMM 
[RUFS2]SPARSE REAL UNSYM 
(GENTAB]GENERATE TWIDDLE 


{CMPACT COMPLEX MATRIX L/U 
(LUF]LU MATRIX 

{ PEEK ] MEMORY 

TO COMPLEX 2-DIMENSIONAL 


TO COMPLEX 2-DIMENSIONAL 
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DOUBLE-PRECISION NUMBER 
DOUBLE-PRECISION SUBTRACT 
DOUBLE-PRECISION SUBTRACT 
EIGENSYSTEM SOLVER 
EIGENSYSTEM SOLVER 
EIGENSYSTEM SOLVER 
EIGENVECTORS 
ELEMENT. .. [ ICAMAX ] INDEX 
ELEMENT IN VECTOR 
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ELEMENT IN VECTOR 
ELEMENT IN VECTOR 
ELEMENT IN VECTOR 
ELEMENT MAGNITUDES 
ELEMENT MAGNITUDES 
ELEMENT SQUARES 
ELEMENT SQUARES 
ELEMENTS 


ELEMENTS... [RMSQV ] ROOT-MEAN-SQUARE 


ELEMENTS 

ELEMENTS INTEGRATION 
EQUAL 
EQUAL...{LVGE]LOGICAL 
EQUAL 

EQUATION SOLVER 
EQUATION SOLVER 
EQUATION SOLVER 
EQUATION SOLVER 
EQUATION, 2 POLES, 2 ZEROS 
EQUIVALENCE 
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EUCLIDEAN NORM 
EVALUATION 

EXCLUSIVE OR 
EXPONENTIAL 
EXPONENTIAL... [CVMEXP ] COMPLEX 
EXPONENTIAL 
EXPONENTIAL AVERAGING 
EXPONENTIAL BASE 19 
FACTOR. ..[CSFR2 ]SPARSE 
FACTOR. ..[CUFR2 ]SPARSE 
FACTOR 

FACTOR 

FACTOR 

FACTOR & SOLVE 

FACTOR & SOLVE 

FACTOR & SOLVE 

FACTOR & SOLVE 

FACTOR TABLE 
PACTORIZATION 
FACTORIZATION CROUT 
FETCH 

FFT...(CFFT2D ]COMPLEX 
FFT...({RFFT2D]REAL 


INDEX 


[GRAD2D ]MAXIMUM 

[ GRD2DB ]MAXIMUM 
[LVGT]LOGICAL VECTOR 
[LVGE]LOGICAL VECTOR 
[VPKR32]VECTOR REAL 

{ VUPR32 | VECTOR 

(CH ]COMPLEX 

[ HTRIBK ] COMPLEX 

[ HTRIDI ]COMPLEX 
{VSINH]VECTOR SINE 
(VTANH ]VECTOR TANGENT 
[COSH]REAL NUMBER 
[SINH]REAL NUMBER 
{TANH]REAL NUMBER 

[ VIMAG ] EXTRACT 
[SCASUM]SUM OF REAL AND 
(MTIMOV]VECTOR MOVE WITH 
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L ae JV GU LUN MOVE WITH 


(TTIMOV]VECTOR MOVE WITH 

{ VINDEX ] VECTOR 
[VSORT]VECTOR SORT WITH 

{ CDOTC ] COMPLEX 

{RKGTF ] R-K-GILL-~THOMPSON 
{ADAMS4]ADAMS VARIABLE STEP 
REAL NUMBER TO NEAREST 

[ VIADD ] VECTOR 

{ VIDIV ] VECTOR 

[ VIFIX]VECTOR 

{ VIMUL ] VECTOR 

{ VINEG 1 VECTOR 
(VPKI32]VECTOR 32-BIT 
{CPOWCI ]COMPLEX TO 
[IPOW]INTEGER TO 
{[RPOWRI]REAL TO 

( VISUB ] VECTOR 

[ VFLOAT ] CONVERT 

(VUSI32 ]VECTOR 32~-BIT SIGNED 
[VISORT]VECTOR SORT 

(RKGIL ]RUNGE-KUTTA-GILL 
SIMPSON'S 1/3 RULE 
[VSUM]VECTOR SUM OF ELEMENTS 
TRAPEZOIDAL RULE 
{SCS1L]SCALAR COS/SIN, TM 
(CFFTI ]COMPLEX FFT WITH 
[NMOLI]NMO LINEAR 
{[NMOQI]NMO QUADRATIC 
(RFFTI]REAL FFT WITH 
{RFTII]REAL FFT WITH QUARTER 
{CMATIN ]COMPLEX MATRIX 
{MATINV ]MATRIX 

{ VRECIP ] VECTOR 

{ PAS2I ]RADIX-2 

{ PAS3I ]JRADIX-3 
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GRADIENT FILTER 

GRADIENT FILTER WITH BOUND 
GREATER THAN 

GREATER THAN OR EQUAL 
HALFWORD PACK 

HALFWORD REAL UNPACK 
HERMITIAN EIGENSYSTEM SOLVER 
HERMITIAN EIGENVECTORS 
HERMITIAN TRIDIAGONALI ZATION 
HYPERBOLIC 

HYPERBOLIC 

HYPERBOLIC 

HYPERBOLIC COSINE 

HYPERBOLIC SINE 

HYPERBOLIC TANGENT 
IMAGINARIES OF COMPLEX VECTOR 
IMAGINARY MAGNITUDES 
INCREMENT MD TO TM 


TNCDEMPEAMT THN TO Mn 
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INCREMENT TM TO TM 

INDEX 

INDICES 

INNER PRODUCT 

INTEG. ORDER 4 

INTEG.ORD 4 

INTEGER... [NINT]ROUND 
INTEGER ADD 

INTEGER DIVIDE 

INTEGER FIX 

INTEGER MULTIPLY 

INTEGER NEGATE 

INTEGER PACK 

INTEGER POWER 

INTEGER POWER 

INTEGER POWER 

INTEGER SUBTRACT 

INTEGER TO FLOATING-POINT 
INTEGER UNPACK 

INTEGER VALUES 
INTEGRATION 

INTEGRATION. ..[VSIMPS ] VECTOR 
INTEGRATION 
INTEGRATION... {VTRAPZ ] VECTOR 
INTERP.ORD 1 
INTERPOLATION 
INTERPOLATION 
INTERPOLATION 
INTERPOLATION 
INTERPOLATION 

INVERSE 

INVERSE 

INVERSE 

INVERSE COMPLEX FFT PASS 
INVERSE COMPLEX FFT PASS 


INDEX - 7 


[ CMMTRC ] COMPLEX 

[ SMPACK ] SPARSE 

{SGESL]REAL GENERAL 

{ LUSN]LU 

{ SGTSL ] TRIDIAGONAL 

[ TRIDIA ] TRIDIAGONAL 

{ SMUPCK ] SPARSE 

[ SMVMUL ] SPARSE 

[ROT3]3D ROTATION 
SYMMETRIC 

{ VMAX ] VECTOR 

[ ISAMAX ] INDEX OF 

{ VMAXMG ] VECTOR 
[MMTMUL]VECTOR MULTIPLY 
[MTMMUL ]VECTOR MULTIPLY 
([MTTMUL]VECTOR MULTIPLY 
[MMTADD]VECTOR ADD 
[MTMADD]VECTOR ADD 
(MT'TADD]VECTOR ADD 
[MMTSUB]VECTOR SUBTRACT 
(MTMSUB]VECTOR SUBTRACT 
{[MTTSUB]VECTOR SUBTRACT 
{[POKE]STORE INTO 

[ TVCLR ] TABLE 
{VLMERG]LOGICAL VECTOR 

[ VMIN ] VECTOR 

{ VMINMG ] VECTOR 

[CVMOV ]COMPLEX VECTOR 
[SVMOV]SPARSE VECTOR 

[ VMOV ] VECTOR 

[ MOVREP ] SUB- IMAGE 
[MTMOV ] VECTOR 

[ TMMOV J VECTOR 
[MTIMOV ] VECTOR 

{ TMIMOV ] VECTOR 

{ TTIMOV ] VECTOR 
[RESNMO]RESIDUAL NORMAL 
[VSMA3]THREE VECTOR SCALAR 
[VSMA4]FOUR VECTOR SCALAR 
(CMVML3 ]COMPLEX 3X3 MATRIX 
(CMVML4 |]COMPLEX 4X4 MATRIX 
[VXCS ] VECTOR 

{ BLKMAN ] BLACKMAN WINDOW 
(CCMMUL ]COMPLEX MATRIX 
[CGMMUL ]COMPLEX GENERAL MATRIX 
[CMMUL ]COMPLEX MATRIX 

[ CMUL ] COMPLEX 

{CRMMUL ]COMPLEX-REAL MATRIX 
AND REAL VECTOR 

[CVMUL ]COMPLEX VECTOR 
[CVSMUL ]COMPLEX VECTOR SCALAR 
TO DOUBLE PRECISION 

TO DOUBLE-PRECISION 

[FMMMV ]FAST MATRIX 
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MATRIX 
MATRIX 
MATRIX 
MATRIX 
MATRIX 
MATRIX 
MATRIX 


MULTIPLY TRACE 

PACK 

SOLVE 

SOLVE CROUT 

SOLVER 

SOLVER 

UNPACK 

MATRIX VECTOR MULTIPLY 
MATRIX, 3-ANGLE 
MATRX...[TRED1]TRIDIAGONALIZE 
MAXIMUM 

MAXIMUM ABSOLUTE VALUE 
MAXIMUM MAGNITUDE 

MD*MD 
MD*TM 
MD*TM 
MD+MD 


9 

2 
88adaa83e 
262252262 


MD-TM 
MEMORY 
MEMORY VECTOR CLEAR 
MERGE 

MINIMUM 

MINIMUM MAGNITUDE 


AND LEVEL REPLACE 

MD TO TM 

TM TO MD 

WITH INCREMENT MD TO TM 
MOVE WITH INCREMENT TM TO MD 
WITH INCREMENT TM TO TM 
MOVEOUT 

MULT AND ADD 

MULT AND ADD 

MULT. 3D VECTORS 

MULT. 4D VECTORS 

MULTIPLIED BY SIN AND COS 
MULTIPLY 
MULTIPLY 
MULTIPLY 
MULTIPLY 
MULTIPLY 
MULTIPLY 
MULTIPLY. 
MULTIPLY 
MULTIPLY 
MULTIPLY... [DMULRR]SINGLE 
MULTIPLY... [DMUL ] DOUBLE 
MULTIPLY 


- - (CRVMUL ] COMPLEX 


INDEX 


[ VNEG ] VECTOR 
[MNAXB ] SUB-MATRIX 
[MNATXB ] SUBMATRIX 

[SCNRM2 ]COMPLEX EUCLIDEAN 

{ SNRM2 ] EUCLIDEAN 

[ RESNMO ] RESIDUAL 
[LVNOT]LOGICAL VECTOR 
[VLNOT]VECTOR LOGICAL 
[LVNE]LOGICAL VECTOR 

[CFFTB ]COMPLEX-TO-COMPLEX FFT 
{RFFTB]REAL-TO-COMPLEX FFT 
(CFFTM ]MIXED-RADIX COMPLEX FFT 
[RFFTM ]MIXED-RADIX REAL FFT 
[AINT]TRUNCATE REAL 

[ATAN ]ARCTANGENT OF REAL 
{[CEXP ]JEXPONENTIAL OF COMPLEX 
[CONJG]CONJUGATE OF COMPLEX 
[CSORT]SQUARE ROOT OF COMPLEX 
[DNEG]NEGATE DOUBLE-PRECISION 
[EXP EXPONENTIAL OF REAL 
[SQRT]SQUARE ROOT OF REAL 

{ ABS ]REAL 

{ CABS ] COMPLEX 

[ ACOS ] REAL 

f ASIN] REAL 

[CCOS ] COMPLEX 

[COS ] REAL 

[RAN]SCALAR RANDOM 

{COSH ] REAL 

[ SINH ] REAL 

{ TANH ] REAL 

[ALOG19 ]REAL 

{ ALOG ] REAL 

{CLOG ] COMPLEX 

{SIGN ]REAL 

{[CSIN ] COMPLEX 

(SIN ]REAL 

[TAN ]REAL 

[NINT]ROUND REAL 

{ANINT]JROUND REAL 

OF RATIO OF REAL 
[VRAND]VECTOR RANDOM 
{FUNL]FUNCTION OF 
[VLOR]VECTOR LOGICAL 

LOGICAL EXCLUSIVE 

VECTOR GREATER THAN 
PREDICTOR 

{AMC1 ]ADAMS-MOULTON CORRECTOR 
PREDICTOR 

[AMC2 ]ADAMS-MOULTON CORRECTOR 
PREDICTOR 

[AMC3 ] ADAMS-MOULTON CORRECTOR 
PREDICTOR 

{AMC4]ADAMS-MOULTON CORRECTOR 
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NEGATE 
NEGATIVE MULTIPLY 

NEGATIVE TRANSPOSE MULTIPLY 
NORM 

NORM 

NORMAL MOVEOUT 

NOT 

NOT 

NOT EQUAL 

NOT IN PLACE 

NOT IN PLACE 

NOT-IN-PLACE 

NOT-IN-PLACE 

NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 
NUMBER 


ABSOLUTE VALUE 

ABSOLUTE VALUE 

ARCCOSINE 

ARCSINE 

COSINE 

COSINE 

GENERATOR 

HYPERBOLIC COSINE 
HYPERBOLIC SINE 
HYPERBOLIC TANGENT 
LOGARITHM 

LOGARITHM 

LOGARITHM 

SIGN TRANSFER 

SINE 

SINE 

TANGENT 

NUMBER TO NEAREST INTEGER 
NUMBER TO NEAREST WHOLE 
NUMBERS... [ATAN2Z ] ARCTANGENT 
NUMBERS 

ONE VARIABLE 

OR 

OR...([VLXOR]VECTOR | 

OR EQUAL...[LVGE]LOGICAL 

ORDER 1...[ABP1]ADAMS-BASHFORTH 
ORDER 1 

ORDER 2...[ABP2 ]ADAMS-BASHFORTH 
ORDER 2 

ORDER 3...[ABP3 ] ADAMS-BASHFORTH 
ORDER 3 

ORDER 4...[ABP4]ADAMS-BASHFORTH 
ORDER 4 
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{ VRAMP ] VECTOR 

[ RAN ] SCALAR 

{ VRAND ] VECTOR 
[POST64]POST BITS TO 
{ATAN2 JARCTANGENT OF 
(MRRUNR ]MIXED-RADIX RFFT 
DOT PRODUCT REAL 

[ SAXPYN ] NESTED 
[SCASUM]SUM OF 

[ SDOTN ] NESTED 

{ RFF'TM ]MIXED-RADIX 
[VPKR32 ] VECTOR 

{ AINT ] TRUNCATE 
[ATAN]ARCTANGENT OF 

[EXP JEXPONENTIAL OF 
[SQRT]SQUARE ROOT OF 

[ NINT ]ROUND 

[ ANINT ] ROUND 

{ATAN2 JARCTANGENT OF RATIO OF 
[CPOWCR ]COMPLEX TO 
[RPOWJREAL TO 
[DDOTRR]DOUBLE DOT PRODUCT 
[RSFS2 ]SPARSE 

(RSFR2 ]SPARSE 

{RSSV2 ]SPARSE 
{VUPR32]VECTOR HALFWORD 
[RUFS2 ]SPARSE 

[RUFR2 ]SPARSE 
[RUSV2 ] SPARSE 

[CRVADD ]COMPLEX AND 
{CRVDIV]COMPLEX AND 
{CRVMUL ]COMPLEX AND 
[CRVSUB]COMPLEX AND 
{[SDOT]DOT PRODUCT OF 
[CVREAL]FORM COMPLEX VECTOR OF 
[ VREAL ] EXTRACT 

[CVRCEP ]COMPLEX VECTOR 
[RRCP ]REAL 

( VRSQRT ] VECTOR 

(RECT ]POLAR TO 
[AMOD]REAL/REAL DIVIDE 
[MOD] INTEGER/INTEGER DIVIDE 
MOVE AND LEVEL 

[ VRVRS ] VECTOR 

{ MRRUNR ]MIXED-RADIX 
[RSQRT]RECIPROCAL SQUARE 
RECIPROCAL SQUARE 
[VSQRT]VECTOR SQUARE 

[ CSQRT ] SQUARE 

[ SORT ] SQUARE 

[CROTG ]COMPLEX GIVENS 
[CSROT]COMPLEX 2-D 
[SROTG]GIVENS PLANE 

{ SROT ] PLANE 
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RAMP . 
RANDOM NUMBER GENERATOR 
RANDOM NUMBERS 

RASTER 

RATIO OF REAL NUMBERS 
RAVEL/UNRAVEL PASS 

REAL... [DDOTRR ] DOUBLE 

REAL A * X+Y 

REAL AND IMAGINARY MAGNITUDES 
REAL DOT PRODUCT 

REAL FFT NOT-IN-PLACE 

REAL HALFWORD PACK 

REAL NUMBER 

REAL NUMBER 

REAL NUMBER 

REAL NUMBER 

REAL NUMBER TO NEAREST INTEGER 
REAL NUMBER TO NEAREST WHOLE 
REAL NUMBERS 

REAL POWER 

REAL POWER 

REAL REAL 

REAL SYMM FACTOR & SOLVE 
REAL SYMMETRIC FACTOR 

REAL SYMMETRIC SOLVE 

REAL UNPACK 

REAL UNSYM FACTOR & SOLVE 
REAL UNSYMMETRIC FACTOR 
REAL UNSYMMETRIC SOLVE 
REAL VECTOR ADD 

REAL VECTOR DIVIDE 

REAL VECTOR MULTIPLY 

REAL VECTOR SUBTRACT 

REAL VECTORS 

REALS 

REALS OF COMPLEX VECTOR 
RECIPROCAL 

RECIPROCAL 

RECIPROCAL SQUARE ROOT 
RECTANGULAR CONVERSION 
REMAINDER 

REMAINDER 

REPLACE. ..{MOVREP ] SUB-IMAGE 
REVERSE ORDERING 
RFFT RAVEL/UNRAVEL PASS 
ROOT 

ROOT. ..{VRSQRT]VECTOR 

ROOT 

ROOT OF COMPLEX NUMBER 
ROOT OF REAL NUMBER 
ROTATION 

ROTATION 

ROTATION 

ROTATION 
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[ VSINH ] VECTOR 

[KSMLV ]K-TH 

COMPLEX SYMM FACTOR AND 
COMPLEX SYMMETRIC 

COMPLEX UNSYM FACTOR AND 
COMPLEX UNSYMMETRIC 

REAL SYMM FACTOR AND 
[RSSV2]SPARSE REAL SYMMETRIC 
REAL UNSYM FACTOR AND 
[RUSV2]SPARSE REAL UNSYMMETRIC 
[SGESL]REAL GENERAL MATRIX 
[LUSN]LU MATRIX 

HERMITIAN EIGENSYSTEM 
MATRIX EQUATION 

MATRIX EQUATION 

[CSOLV ]COMPLEX SYSTEM 
SYMMETRIC EIGENSYSTEM 
[RS]REAL SYMMETRIC EIGENSYSTEM 
[SGTSL]TRIDIAGONAL MATRIX 
[SITSOL]SPARSE ITERATIVE 
FORMAT EQUATION 
[SOLVEQ]LINEAR EQUATION 
[TRIDIA]TRIDIAGONAL MATRIX 
[ VASORT ] VECTOR 

[ VISORT ] VECTOR 

[ VSORT ] VECTOR 

[VSQ] VECTOR 

[VSSQ]VECTOR SIGNED 

[RSQRT ]RECIPROCAL 
[VRSQRT]VECTOR RECIPROCAL 
[ VSORT ] VECTOR 

VECTOR MAGNITUDE 
[MEASQV]MEAN OF VECTOR ELEMENT 
[SVESQ]SUM OF VECTOR ELEMENT 
[SVS]SUM OF VECTOR SIGNED 
[ADAMS4]ADAMS VARIABLE 

[ CMTRAC ] COMPLEX 

[ CMTRAN ] COMPLEX 

AND REAL VECTOR 
[CVSUB]COMPLEX VECTOR 

TO DOUBLE-PRECISION 

TO DOUBLE-PRECISION 
(VISUB]VECTOR INTEGER 
MULTIPLY, MULTIPLY, AND 
[VMSB]VECTOR MULTIPLY AND 
SCALAR MULTIPLY AND 
[VSUB ] VECTOR 

[ VSBM] VECTOR 

[ VSBSM ] VECTOR 

[ MMTSUB ] VECTOR 

[ MTMSUB ] VECTOR 

[ MT'TSUB ] VECTOR 

[ TMMSUB ] VECTOR 

[ TIMTSUB ] VECTOR 


FPS 869-7482-991C 


SINE HYPERBOLIC 
SMALLEST ELEMENT IN VECTOR 
SOLVE. ..[CSFS2 ]SPARSE 
SOLVE. ..[CSSV2 ]SPARSE 
SOLVE. ..[CUFS2 ]SPARSE 
SOLVE... [CUSV2 ]SPARSE 
SOLVE. ..[{RSFS2 ]SPARSE 
SOLVE 

SOLVE... [RUFS2 ]SPARSE 
SOLVE 

SOLVE 

SOLVE CROUT 
SOLVER... [CH] COMPLEX 
SOLVER. .. { CMSOLV ] COMPLEX 
SOLVER... {CSOLVQ ] COMPLEX 
SOLVER 

SOLVER... [EIGRS ]REAL 
SOLVER 

SOLVER 

SOLVER — 

SOLVER... [{SKYSOL]SKYLINE 
SOLVER 

SOLVER 

SORT ALGEBRAIC VALUES 
SORT INTEGER VALUES 

SORT WITH INDICES 

SQUARE 

SQUARE 

SQUARE ROOT 

SQUARE ROOT 

SQUARE ROOT 
SQUARED... [CVMAGS ] COMPLEX 
SQUARES 

SQUARES 

SQUARES 

STEP INTEG.ORD 4 
SUB~-MATRIX TRACE 
SUB-MATRIX TRANSPOSE 
SUBTRACT. ..[CRVSUB ]COMPLEX 
SUBTRACT 
SUBTRACT. . . { DSUBRR ] SINGLE 
SUBTRACT. ..[DSUB]DOUBLE 


SUBTRACT 

SUBTRACT. ..[VMMSB] VECTOR 
SUBTRACT 
SUBTRACT. ..[ VSMSB] VECTOR 
SUBTRACT 


SUBTRACT AND MULTIPLY 
SUBTRACT AND SCALAR MULTIPLY 
SUBTRACT MD-MD TO TM 
SUBTRACT MD-TM TO MD 
SUBTRACT MD-TM TO TM 
SUBTRACT TM-MD TO MD 
SUBTRACT TM-MD TO TM 
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(TTTADD]VECTOR ADD 
[TMMSUB]VECTOR SUBTRACT 
[TMTSUB]VECTOR SUBTRACT 
[TTMSUB]VECTOR SUBTRACT 
[TTTSUB]VECTOR SUBTRACT 

MATRIX MULTIPLY 

[CMTRAC ]COMPLEX SUB-MATRIX 
(ISIGN] INTEGER SIGN 
{SIGN]REAL NUMBER SIGN 
[CTRN2]2-D COORDINATE 
[RDFT]REAL DISCRETE FOURIER 
COORDINATE 

{HLBRT ]HILBERT 

(CMTRAN ]COMPLEX SUB-MATRIX 
[MAXBT ]MATRIX A TIMES B 
[MTRANS ]MATRIX 

[ STMM ] SUBMATRIX 

[MATXBT ]SUBMATRIX TRANSPOSE 
{[MNATXB]SUBMATRIX NEGATIVE 
[MATXBT ] SUBMATRIX 

[ VTRAPZ ] VECTOR 

[V@1L]VECTOR ZERO 

[ IMTQL1 ] DIAGONALIZE 

{ IMTQL2 ]DIAGONALIZE A 
[HTRIDI ]COMPLEX HERMITIAN 

[ VAINT ] VECTOR 

[ VFRAC ] VECTOR 

[ GENTAB ]GENERATE 

[ SMUPCK]SPARSE MATRIX 
{SVUPCK]SPARSE VECTOR 
(VUP16]VECTOR 16-BIT BYTE 
(VUP32]VECTOR 32-BIT BYTE 
(VUP8]VECTOR 8-BIT BYTE 
(VUPR32]VECTOR HALFWORD REAL 
16-BIT SIGNED BYTE 

32-BIT SIGNED BYTE 

8-BIT SIGNED BYTE 

32-BIT SIGNED INTEGER 
[VUUI32]VECTOR 32-BIT UNSIGNED 
[VUUI32]VECTOR 32-BIT 
(CUFS2]SPARSE COMPLEX 
[RUFS2]SPARSE REAL 
[CUFR2]SPARSE COMPLEX 
[RUFR2]SPARSE REAL 
[CUSV2]SPARSE COMPLEX 
[RUSV2]SPARSE REAL 
[CCEPS ] PHASE 
{SHPHU]SCHAFER'S PHASE 
[PKVAL]PEAK AND 
[VASORT]VECTOR SORT ALGEBRAIC 
[VISORT]VECTOR SORT INTEGER 
(FUNL]FUNCTION OF ONE 

{ ADAMS4 ] ADAMS 

[VARNMO]NMO WITH 


FPS 869-7482-S91C 


q 


4 KEY WORD INDEX 


TM-T TO TM : 
TRACE... [CMMTRC ] COMPLEX 
TRACE 

TRANSFER 

TRANSFER 

TRANSFORM 

TRANSFORM 
TRANSFORMATION. ..[{CTRN3 ] 
TRANSFORMER 

TRANSPOSE 

TRANSPOSE 

TRANSPOSE 

TRANSPOSE & MULTIPLY 
TRANSPOSE MULTIPLY 
TRANSPOSE MULTIPLY 
TRANSPOSE TRANSPOSE MULTIPLY 
TRAPEZOIDAL RULE INTEGRATION 
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VALUES 

VALUES 

VARIABLE 
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[ SSWAP ] INTERCHANGES 
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[ BLKMAN ] BLACKMAN 

[ HAMM ] HAMMING 
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VECTORS ; 

VECTORS. ..[{SN2]SQUARED 
VECTORS 

VELOCITY 

VELOCITY 
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WINDOW MULTIPLY 

WINDOW MULTIPLY 
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WORKSPACE 

X+yY 


++ + 


Y 
XY 
XY 


Ww KK OM OM OM 


ZERO TRENDS 
ZEROS... [DEQ22 ] DIFFERENCE 
ZEROS 


INDEX 


Please detach cards along perforations. 


Your comments will help us improve the quality and usefulness of our publications. Please fill 
out and return this form. (The mailing address is on the back.) 


Title of document: 
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Firm: Department: 
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City: State: Zip Code: 
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I used this manual. . . I found this material. . . 

Yes No 
QI) as an introduction to the subject accurate WY 
CJ as an aid for advanced training complete Oo 
QJ to instruct a class written clearly O QO 
OC to learn operating procedures well illustrated oO QO 
CY as a reference manual well indexed UO O 
OD other 


Please indicate below, listing the pages, any errors you found in the manual. Also indicate if 
you would have liked more information about a certain subject. 
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